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Abstract 

The classical Euler's problem on stationary configurations of elastic 
rod with fixed cndpoints and tangents at the endpoints is considered as a 
left-invariant optimal control problem on the group of motions of a two- 
dimensional plane E(2). 

The attainable set is described, existence and boundedness of opti- 
mal controls are proved. Extremals are parametrized by Jacobi's elliptic 
functions of natural coordinates induced by the flow of the mathematical 
pendulum on fibers of the cotangent bundle of E(2). 

The group of discrete symmetries of Euler's problem generated by re- 
fiections in the phase space of the pendulum is studied. The corresponding 
Maxwell points are completely described via the study of fixed points of 
this group. As a consequence, an upper bound on cut points in Euler's 
problem is obtained. 
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1 Introduction 



In 1744 Leonhard Euler considered the following problem on stationary configu- 
rations of elastic rod. Given a rod in the plane with fixed endpoints and tangents 
at the endpoints, one should determine possible profiles of the rod under the 
given boundary conditions. Euler obtained ODEs for stationary configurations 
of the elastic rod and described their possible qualitative types. These configu- 
rations are called Euler elasticae. 

An Euler elastica is a critical point of the functional of elastic energy on the 
space of smooth planar curves that satisfy the boundary conditions specified. 
In this paper we address the issue of optimality of an elastica: whether a critical 
point is a minimum of the energy functional? That is, which elasticae provide 
the minimum of the energy functional among all curves satisfying the boundary 
conditions (the global optimality), or the minimum compared with sufficiently 
close curves satisfying the boundary conditions (the local optimality). These 
questions remained open despite their obvious importance. 

For the elasticity theory, the problem of local optimality is essential since 
it corresponds to stability of Euler elasticae under small perturbations that 
preserve the boundary conditions. In the calculus of variations and optimal 
control, the point where an extremal trajectory loses its local optimality is called 
a conjugaie point. We will give an exact description of conjugate points in the 
problem on Euler elasticae, which were previously known only numerically. 

From the mathematical point of view, the problem of global optimality is 
fundamental. We will study cut points in Euler's elastic problem — the points 
where elasticae lose their global optimality. 

This is the first of two planned works on Euler's elastic problem. The aim 
of this work is to give a complete description of Maxwell points, i.e., points 
where distinct extremal trajectories with the same value of the cost functional 
meet one another. Such points provide an upper bound on cut points: an ex- 
tremal trajectory cannot be globally optimal after a Maxwell point. In the 
second work [34] we prove that conjugate points in Euler's elastic problem are 
bounded by Maxwell points. Moreover, we pursue the study of the global opti- 
mal problem: we describe the global diffeomorphic properties of the exponential 
mapping. 

This paper is organized as follows. In Sec. 2 we review the history of the 

problem on elasticae. In Sec. 3 we state Euler's problem as a left-invariant 
optimal control problem on the group of motions of a two-dimensional plane E(2) 
and discuss the continuous symmetries of the problem. In Sec. 4 we describe the 
attainable set of the control system in question. In Sec. 5 we prove existence 
and boundedness of optimal controls in Euler's problem. In Sec. 6 we apply 
Pontryagin Maximum Principle to the problem, describe abnormal extremals, 
and derive the Hamiltonian system for normal extremals. 

Due to the left-invariant property of the problem, the normal Hamiltonian 
system of PMP becomes triangular after an appropriate choice of parametriza- 
tion of fibers of the cotangent bundle of E(2): the vertical subsystem is in- 
dependent on the horizontal coordinates. Moreover, this vertical subsystem is 
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essentially the equation of the mathematical pendulum. For the detailed subse- 
quent analysis of the extremals, it is crucial to choose convenient coordinates. 
In Sec. 7 we construct such natural coordinates in the fiber of the cotangent 
bundle over the initial point. First we consider the "angle- action" coordinates 
in the phase cylinder of the standard pendulum, and then continue them to the 
whole fiber via continuous symmetries of the problem. One of the coordinates is 
the time of motion of the pendulum, and the other two are integrals of motion 
of the pendulum. In Sec. 8 we apply the elliptic coordinates thus constructed 
for integration of the normal Hamiltonian system. In particular, we recover the 
different classes of elasticae discovered by Leonhard Euler. 

The flow of the pendulum plays the key role not only in the paxametrization 
of extremal trajectories, but also in the study of their optimality. In Sec. 9 we 
describe the discrete symmetries of Euler's problem generated by reflections in 
the phase cylinder of the standard pendulum. Further, we study the action of 
the group of reflections in the preimage and image of the exponential mapping 
of the problem. 

In Sec. 10 we consider Maxwell points of Euler's problem. The Maxwell 
strata corresponding to reflections are described by certain equations in elliptic 
functions. In Sec. 11 we study solvability of these equations, give sharp estimates 
of their roots, and describe their mutual disposition via the analysis of the 
elliptic functions involved. 

A complete description of the Maxwell strata obtained is important both for 
global and for local optimality of extremal trajectories. In Sec. 12 we derive 
an upper bound on the cut time in Euler's problem due to the fact that such a 
trajectory cannot be globally optimal after a Maxwell point. In our subsequent 
work [34] we will show that conjugate points in Euler's problem are bounded by 
Maxwell points and give a complete solution to the problem of local optimality 
of extremal trajectories. 

In Sec. 13 we collect the definitions and properties of Jacobi's elliptic func- 
tions essential for this work. 

We used the system "Mathematica" [39] to carry out complicated calcula- 
tions and to produce illustrations in this paper. 

Acknowledgment. The author wishes to thank Professor A.A.Agrachev for bring- 
ing the pearl of Euler's problem to author's attention, and for numerous fruitful 
discussions of this problem. 

2 History of Euler's elastic problem 

In addition to the original works of the scholaxs who contributed to the the- 
ory, in this section we follow also the sources on history of the subject by 
C.Truesdell [37], A.E.H.Love [24], and S.Timoshenko [36]. 

In 1691 James Bernoulli considered the problem on the form of a uniform 
planar elastic bar bent by the external force F. His hypothesis was that the 

B 

bending moment of the rod is equal to — , where B is the "flexural rigidity". 
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and R is the radius of curvature of bent bar. For elastic bar of unit excursion 
built vertically into a horizontal wall and bent by a load sufficient to make its 
top horizontal (rectangular elastica, see Fig. 1), James Bernoulli obtained the 
ODEs 

x"^ dx , dx ^ , 

(where (x,y) is the clastic bar, and s is its length parameter), integrated them 
in series and calculated precise upper and lower bounds for their value at the 
endpoint x = 1, see [8]. 




• x 

1 



Figure 1: James 
Bernoulli's rectan- 
gular elastica 



In 1742 Daniel Bernoulli in his letter [7] to Leonhard Euler wrote that the 
elastic energy of bent rod is proportional to the magnitude 

and suggested to find the elastic curves from the variational principle E min. 

At that time Euler was writing his treatise on the calculus of variations 
"Mcthodus invcnicndi ..." [13] published in 1744, and he adjoined to his book an 
appendix "De curvis elasticis" , where he applied the newly developed techniques 
to the problem on elasticae. Euler considered a thin homogeneous elastic plate, 
rectilinear in the natural (unstressed) state. For the profile of the plate, Euler 
stated the following problem: 

"That among all curves of the same length which not only 

pass through the points A and B, hut are also tangent 

to given straight lines at these points, that curve he determined (2.1) 

/•^ ds 

in which the value of / ^ he a minimum. " 

J A 
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Eulcr wrote down the ODE known now as Euler-Lagrange equation for the 
corresponding problem of calculus of variations and reduced it to the equations: 



dy = 



(a + (3x + 72;^) dx 
sja^ - (a + + 7a;2)2 



ds = 



a? dx 



-{a + I3x + 7x2)2 ' 



(2.2) 



a /3 

where — and 7 are real parameters expressible in terms of B, the load 
? . 

of the elastic rod, and its length. Euler studied the quadrature determined 
by the first of equations (2.2). In the modern terminology, he investigated 
the qualitative behavior of the elliptic functions that parametrize the elastic 
curves via the qualitative analysis of the determining ODEs. Euler described 
all possible types of elasticae and indicated the values of parameters for which 
these types are realized (see a copy of Euler's original sketches at Fig. 2). 




Figure 2: Euler's sketches 



Euler divided all elastic curves into nine classes, they are plotted respectively 

as follows: 

1. straight line. Fig. 16, 

2. Fig. 17, 
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3. rectangular elastica, Fig. 18, 



4. Fig. 19, 

5. periodic elastica in the form of figure 8, Fig. 20, 

6. Fig. 21, 

7. elastica with one loop, Fig. 22, 

8. Fig. 23, 

9. circle, Fig. 24. 

Following the tradition introduced by A. E.H.Love [24], the elastic curves 
with inflection points (classes 2-6) are called inflectional, the elastica of class 
7 is called critical, and elasticae without inflection points of class 8 are called 
non-inflectional. 

Further, Euler established the magnitude of the force applied to the elas- 
tic plate that results in eac;li type of {Jastic:ae. He indicated the experimental 
method for evaluation of the flexural rigidity of the elastic plate by its form un- 
der bending. Finally, he studied the problem of stability of a column modeled by 
the loaded rod whose lower end is constrained to remain vertical, by presenting 
it as an elastica of the class 2 close to the straight line (thus a sinusoid). 

After the work of Leonhard Euler, the elastic curves are called Euler elasti- 
cae. 

The first explicit parametrization of Euler elasticae was performed by L.Saal- 
chiitz in 1880 [29]. 

In 1906 the future Nobel prize-winner Max Born defended a Ph.D. thesis 
called "Stability of elastic lines in the plane and the space" [10]. Born considered 
the problem on elasticae as a problem of calculus of variations and derived from 
Euler-Lagrange equation that its solutions {x{t),y{t)) satisfy the ODEs of the 
form: 

X = cos 9, y = sin 9, 

Ae -\- Rsin{0 - j) = 0, A, i?, 7 = const, (2.3) 

thus the angle 9 determining the slope of elasticae satisfies the equation of the 
mathematical pendulum (2.3). 

Further, Born studied stability of elasticae with fixed endpoints and fixed 
tangents at the endpoints. Born proved that an elastic; arc; without inflection 
points is stable (in this case the angle 6 is monotone, thus it can be taken 
as a parameter along elastica; Born showed that the second variation of the 

functional of elastic energy E = ^ J 6^ dt is positive). In the general case. 

Born wrote down the Jacobian that vanishes at conjugate points. Since the 
functions entering this Jacobian were too complicated. Born restricted himself 
to numerical investigation of conjugate points. He was the first to plot elasticae 
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numerically and check the theory against experiments on elastic rods, see the 
photos from Born's thesis at Fig. 3. Moreover, Born studied stability of Euler 
elasticae with various other boundary conditions, and obtained some results for 
elastic curves in M.^. 




Figure 3: Max Born's experiments 

In 1986 A.Arthur and G.R.Walsh [6] and, independently, in 1993 V.Jurdje- 
vic [19] discovered that Euler elasticae appear in the ball-plate problem stated 
as follows. Consider a ball rolling on a horizontal plane without slipping or 
twisting. The problem is to roll the ball from an initial contact configuration 
(defined by contact point of the ball with the plane, and orientation of the ball 
in the 3-space) to a terminal contact configuration, so that the curve traced by 
the contact point in the plane was the shortest possible. Arthur and Walsh, 
and Jurdjevic showed that such optimal curves are Euler elasticae. Moreover, 
Jurdjevic also extensively studied the elastic problem in M'^, its analogs in the 
sphere S^, and in the Lorentz space [20,21]. 

In 1993 R.Brockett and L.Dai [11] discovered that Euler elasticae are pro- 
jections of optimal trajectories in the nilpotent sub-Riemannian problem with 
the growth vector (2,3,5) known also as generalized Dido problem [30-33]. 

Elasticae were considered in approximation theory as nonlinear splines [9, 
14,17,18,23], in computer vision as a maximum likelihood reconstruction of 
occluded edges [27], their 3-dimensional analogues are used in the modeling of 
DNA minicircles [25, 26] etc. 

Euler elasticae and their various generalizations play an important role in 
modern mathematics, mechanics, and their applications. Although, the initial 
variational problem as it was stated by Euler (2.1) is far from complete solution: 
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neither local nor global optimality of Euler elasticae is studied. This is the first of 
two planned works that will give a complete description of local optimality, and 
present an essential progress in the study of the global optimality of elasticae. In 
this paper wc give an upper bound on the cut points along Euler elasticae, i.e., 
points where they lose their global optimality. In the next work [34] we obtain 
a complete characterization of conjugate points, i.e., points where elasticae lose 
their local optimality. 

We would like to complete this historical introduction by two phrases of 
S.Antman [4]. On the one hand, "Fortunately Euler left some unsolved issues 
for his successors," but on the other hand, "There is unfortunately a voluminous 
and growing literature devoted to doing poorly what Euler did well." With the 
hope to contribute to the first tradition rather than to the second one, we start 
this work. 



3 Problem statement 
3.1 Optimal control problem 

First wc state the elastic problem mathematically. Let a homogeneous elastic 
rod in the two-dimensional Euclidean plane have a fixed length I > 0. Take 
any points oq, ai G and arbitrary unit tangent vectors at these points 
Vi G TaM'^, \vi\ = 1, i = 0,1. The problem consists in finding the profile of a 
rod 7 : [0,ti] M^, starting at the point qq and coming to the point ai with 
the corresponding tangent vectors Vq and vi: 

7(0) -flo, 7(ii)-ai, (3.1) 
7(0)=t;o, i{ti)=vi, (3.2) 

with the minimal elastic energy. The curve j{t) is assumed absolutely contin- 
uous with Lebesgue square-integrable curvature k{t). We suppose that j{t) is 
arc-length parametrized, i.e., \'y{t)\ = 1, so the time of motion along the curve 
7 coincides with its length: 

ti = I. (3.3) 
The elastic energy of the rod is measured by the integral 




We choose Cartesian coordinates {x, y) in the two-dimensional plane K^. Let 
the required curve be parameterized as j(t) = {x{t),y(t)), t G [0,ii], and let 
its endpoints have coordinates = {xi,yi), i = 0, 1. Denote by 9 the angle 
between the tangent vector to the curve 7 and the positive direction of the 
axis X. Further, let the tangent vectors at the endpoints of 7 have coordinates 
Vi = (cos ^i, sin ^j), i = 0,1, see Fig. 4. 
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y 




X 



Figure 4: Statement of Euler's prob- 
lem 



Then the required curve 7(f) = {x{t),y{t)) is determined by a trajectory of 
the following control system: 

X = cos e, (3.4) 

y = sin6', (3.5) 

e = u, (3.6) 

q = (x, y, 6) e M = M^. ,^ x , u e E, (3.7) 

g(0) =qo = {xQ,yo,0Q), q{ti) = qi = {xi,yi,0i), h fixed. (3.8) 

For an arc-length parametrized curve, the curvature is, up to sign, equal to 
the angular velocity: k"^ = 9^ = u^, whence we obtain the cost functional 

1 Z"*^ 

J = - u^{t)dt^ min . (3.9) 

2 Jo 

We study the optimal control problem (3.4)-(3.9). Following V.Jurdjevic [21], 
this problem is called Euler's elastic problem. Admissible controls are u{t) e 
L2[0,ti], and admissible trajectories are absolutely continuous curves q{t) e 
AC{[QM]\M). 

In vector notations, the problem reads as follows: 

q = Xi{q) + uX2{q), q€M = R^xS\ u € R, (S) 
q{0) = qo, qiti) = qi, t\ fixed, 



1 



J = TZ I u^dt min, 
u e L2[0,ti], 

where the vector fields in the right-hand side of system S are: 

d d d 

Xi = cos6'— +sin6'— , = 

ox oy of) 
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Notice the multiplication table in the Lie algebra of vector fields generated by 
Xi, X2: 

[Xi, X2] = X3 = sin 0— - cos 9—, (3.10) 

ox ay 

[X2,X3]=Xi, [Xi,X2]=0. (3.11) 

3.2 Left-invariant problem 

on the group of motions of a plane 

Euler's elastic problem has obvious symmetries — parallel translations and ro- 
tations of the two-dimensional plane M?. Thus it can naturally be stated as an 
invariant problem on the group of proper motions of the two-dimensional plane 

cos 9 — sin c 
E(2)=.{| sin6» cos6' y | | (x, y) e M^ 6» e 5^ 
1 

Indeed, the state space of the problem M = x Sg is parametrized by 
matrices of the form 



eE(2), 

and dynamics (3.4)-(3.6) is left-invariant on the Lie group E(2): 

^ I — a. 1 / — ttsin^ — ucos^ cos^ 

UC0&9 —usm9 sin 6* 






cos 


t) — sm 


sin 1 


9 cos 6 








e 


— sin^ 


9 


cos 9 











The Lie algebra of the Lie group E(2) has the form 

e(2) = span(£'2i - -B12, Eis, E23), 

where Eij denotes the 3x3 matrix with the only identity entry in the i-th row 
and j-th column, and zeros elsewhere. In the basis 

ei = Eis, 62 = E21 — E12, 63 = —E23, 

the multiplication table in the Lie algebra e(2) takes the form 

[61,62] = 63, [62,63] = 61, [61,63] =0. 
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Then Eulcr's clastic problem becomes the following left-invariant problem 
on the Lie group E(2): 

q = Xi{q)+uX2{q), 9 G E(2), uGR. 
a(0) = qo: q{ti) = 91, h fixed, 

J = - / u^dt min, 
2 Jo 

where 

X,{q) = qe„ i=l,2, g € E(2), 

are basis left-invariant vector fields on E(2) (here gcj denotes the product of 
3x3 matrices). 

3.3 Continuous symmetries 

and normalization of conditions of the problem 

Left translations on the Lie group E(2) are symmetries of Euler's elastic problem. 
By virtue of these symmetries, we can assume that initial point of trajectories 
is the identity element of the group 




I.e., 

r?o = (xo, 2/0, 0o) = (0,0,0). (3.12) 

In other words, parallel translations in the plane ^x,y ^^^^^ the initial point of 
the elastic rod 7 to the origin (0, 0) S and rotations of this plane combine 

the initial tangent vector 7(0) with the positive direction of the axis x. 

Moreover, one can easily see one more continuous family of symmetries of the 
problem — dilations in the plane ^x,y Consider the following one-parameter 
group of transformations of variables of the problem: 

{x, y, 6,t, u, ti, J) I— > {x, y, u, t\,J) = {e^x, e^y, 9, e^t, e~^u, e^t\,e~'^ J). 

(3.13) 

One immediately checks that Euler's problem is preserved by this group of 
transformations. Thus, choosing s = — Inii, we can assume that ti = 1. In 
other words, we obtain an elastic rod of unit length by virtue of dilations in the 
plane R^. ^. 

In the sequel we usually fix the initial point qo as in (3.12). Although, the 
terminal time ti will remain a parameter, not necessarily equal to 1. 

4 Attainable set 

Consider a smooth control system of the form 

9 = /(«,«), 9eM, ueU. (4.1) 
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Let u = u(t) be an admissible control, and let qo G M. Denote by q(t: u, qo) 
the trajectory of the system corresponding to the control u{t) and satisfying the 
initial condition q{0;u,qo) = qo- Attainable set of system (4.1) from the point 
go for time is defined as follows: 

Aqg{ti) = {q{ti;u,qo) \ u = u{t) admissible control, t G [0,ti]}. 

Moreover, one can consider the attainable set for time not greater than ti: 

0<t<ti 

and the attainable set for arbitrary nonnegative time: 

"^90 ~ U "^90 (^)- 

0<t<oo 

The orbit of the system (4.1) is defined as 

^90 = {e""^" o . . . o e^^^Hao) InGR, fi = f{;Ui), meU, NeN}, 

where e'^'^' is the flow of the vector field See [2,21] for basic properties of 
attainable sets and orbits. 

In this section we describe the orbit and attainable sets for Euler's elastic 
problem. 

Multiplication rules (3.10), (3.11) imply that control system S is full-rank; 

Lieq(Xi,X2) = span(Xi(g),X2(g),X3(g)) = T,M V g e M. 

By the Orbit Theorem of Nagano-Sussmann [2,21], the whole state space is a 
single orbit: 

Og„=M yqoGM. 
Moreover, the system is completely controllable: 

Aq,=M yqoGM. 

This can be shown either by appljdng a general controllability condition for 
control-affine systems with recurrent drift (Th. 5 in Sec. 4 [21]), or via con- 
trollability test for left-invariant systems on semi-direct products of Lie groups 
(Th. 10 in Sec.6 [21]). 

On the other hand, it is obvious that system E is not completely controllable 
on a compact time segment [0,fi]: 

All ^ M 

in view of the bound {x{t) — xq)^ + {y{t) — yo)^ < t\, the distance between the 
endpoints of the elastic rod should not exceed the length of the rod. We have 
the following description of the exact-time attainable sets for Euler's problem. 
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Theorem 4.1. Let go = {xo,yo,Oo) S M = x 5^ and ti > 0. Then the 

attainable set of system S is 

Aoih) = {{x, y,0)eM\{x- xo)2 + {y- yaf < tj 

or {x,y,0) = {xo + ticos9o,yo + tism0o,9o)}. 

Proof. In view of the continuous symmetries of the problem (see Subsec. 3.3), 
it suffices to prove this theorem in the case qo = M = (0,0,0), ti = 1, so we 
show that 

^ = ^id(l) = {{x,y,9)GM\x'' + y^ < 1 or {x,y,9) = (1,0,0)}. 

(1) It is easy to see that 

x^ + y'^>l ^ q={x,y,9)^A. 

Indeed, the curve 7(f) = {x{t),y{t)) has unit velocity, thus 

x\t)+y'it) = m-'y{0)\'<t\ (4.2) 

So points {x, y) at the distance greater then 1 from the origin 7(0) = (0, 0) are 
not attainable from the origin for time 1. 

(2) Let + y2 = 1. We show that 

q={x,y,9)GA <^ g= (1,0,0). 

It is obvious that the point q = (1, 0, 0) is attainable from the point qo = (0, 0, 0) 

for time 1 via the control u{t) = 0, t E [0, 1]. 

Conversely, let q = {x,y,9) G A. Consider the function f{t) = x^{t) + 
y'^{t) e W2,2, where q{t) = {x{t),y{t), 6{t)), t G [0, 1], is a trajectory of system S 
connecting the points go and q. We prove that f{t) = t^ , t G [0,1]. It was 
shown in (4.2) that f{t) < t^, t & [0,1]. The function f{t) takes the same 
values as at the endpoints of the segment [0,1]. So if /(t) ^ t^, t e [0, 1], 
then /'(to) > {t^y\t=ta = 2^0 at some point Iq G [0,1]. But this inequality is 
impossible in view of the chain 

\f'it)\ = 2\xx + yy\ < 2^/x^ + y^ = 2^/J{i) <2t, t€ [0, 1]. 

Hence f(t) = t'^, t G [0,1], and the preceding inequalities turn into equalities. 
Then {x{t),y{t)) = a{t){x{t),y{t)), a{t) > 0, whence it follows that ^ = and 
g= (1,0,0). 

(3) Finally, we show that for any angle 9 G 

x^ + y^<l q={x,y,e)GA. 

First we mention some simple trajectories of system S. In the case u = 
we obtain a straight line {x{t),y{t)), and in the case u = C ^ the curve 
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{x{t),y{t)) is a circle of radius — — . Notice that the time of complete revolution 

ici 

27r 

along such circle is j- — ^ as C — > oo. 

Now we construct a trajectory of system E connecting the initial point qo = 
(0, 0, 0) with the terminal one q = {x, y, 6), x"^ + y"^ < 1, for time 1. 

Assume first that {x, y) ^ (0, 0). In the plane y, construct a circle of small 
radius starting at the point (0,0) with the tangent vector (1,0), and a circle of 
small radius starting at the point {x^y) with the tangent vector (cos 0, sin 0). It 
is obvious that there exists a straight line segment in the plane y tangent 
at its initial point to the first circle, and tangent at its terminal point to the 
second circle, in such a way that direction of motion along the circles and the 
segment was consistent, see Fig. 5. 




M = C3 



Figure 5: Steering go to q 



In such a way we obtain a trajectory of system S corresponding to a piecewise 
constant control u(-) taking some values Ci ^ 0, 0, C2 7^ 0; this trajectory 
projects to the plane y to the concatenation of a circle arc, a line segment, 

and a circle arc. Choosing circles of sufficiently small radii j^j~y J^j~y 

obtain the total time of motion along this trajectory t < 1. In order to have a 
trajectory with the same endpoints at the time segment i £ [0, 1], it is enough 

1 1 - T 

to add a circle of radius — = before the first circle, see Fig. 5. The 

63 2n 

trajectory constructed steers the point qa = (0,0,0) to the point q = {x,y,9) 
for time 1. 

If (x, y) = (0,0), then we move from the point (0,0) along a short segment 
to a point (e, 0), and repeat the preceding argument. 

Now the statement of Th. 4.1 follows from the statements (l)-(3) proved 
above. □ 

The following properties of attainable sets of system E follow immediately 
from Th. 4.1. 

Corollary 4.1. Let go be an arbitrary point of M . Then: 
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(1) Agoih) C ^qo(*2) for any < ii < t2. 

(2) Al^=Ag,it) foranyt>0. 

(3) qo e int A^^ for any t>0. 

Item (3) means that system E is small-time locally controllable. Although, 

the restriction of E to a small neighborhood of a point qq E M is not controllable 
since some points in the neighborhood of go are reachable from go by trajectories 
of E far from go- 

Topologically. the attainable set Aq„{t) is an open solid torus united with 
a single point at its boundary. In particular, the attainable set is neither open 
nor closed. 

In the sequel we study Euler's problem under the natural condition 

gie^,„(ii). (4.3) 

5 Existence and regularity of optimal solutions 

We apply known results of optimal control theory in order to show that in 
Euler's elastic problem optimal controls exist and are essentially bounded. 

5.1 Embedding the problem into 

The state space and attainable sets of Euler's problem have nontrivial topology, 
and we start from embedding the problem into Euclidean space. By Th. 4.1, 
the attainable set A = AgQ{l), go = (0,0,0), is contained in the set 

M = cM = {{x, y,e)&M\x^ + y'^ < 1}. 

Moreover, by item (2) of Corollary 4.1, any trajectory of system E starting at 
go does not leave the set M at the time segment f G [0, 1]. So this set can be 
viewed as a new state space of the problem. The set M is embedded into the 
Euclidean space R^^^^^^ by the diffeomorphism 

(j) • M ^ ]R3 

^{x,y,9) = {xi,X2,X3) = {{2 + x)cos0,{2 + x)sia0,y). (5.1) 

The image 

$(M) = {{xuX2,X3) eR^\{2-pf + xl<l}, p= ^Jxl + xl, 
is the closed solid torus. 
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In the coordinates {xi,X2,xs), Euler's problem reads as follows: 

x"^ 

Xi = 3 2 2 ~ "^2, (5.2) 

a=2 = -2— -2 (5.3) 

ti/-]^ ~l~ "-"2 

X = (xi, X2, xz) e $(M), M G K, (5.5) 
x(0)=a;° = (2,0,0), a;(l) = = (a;l,x^,xj), (5.6) 

1 /"^ 

J = - / w^dt ^ min, (5.7) 

2 io 

u(-) e L2[0, 1], a;(-) e AC[0, 1]. (5.8) 
5.2 Existence of optimal controls 

First we cite an appropriate general existence result for control-affine systems 
from Sec. 11.4.C of the textbook by L.Cesari [12]. Consider optimal control 
problem of the form: 

m 

x= j{t,x)+^^Uig,(t,x), a;eXcK", u = (mi, ...,«„)€ M"*, (5.9) 



1=1 



/■ti 

— \ f^^i ^' ^) ~^ "csmi^ 

JQ 



(5.10) 



x{-)eAC{[0,h],X), u(-) eL2([0,ii],M"), (5.11) 
a;(0) = a;°, x{ti) = x\ h fixed. (5.12) 

For such a problem, a general existence theorem is formulated as follows. 

Theorem 5.1 (Th. 11.4.VI [12]). Assume that the following conditions hold: 

(C) the set X is closed, and the function fo is continuous on [0, ti]x X x R™, 

(Li) there is a real-valued function ^{t) >0,t€ [0,ti], ip G Li[0,ii], such that 
fo{t,x,u) > -il>{t) for {t,x,u) e [0,ii] x X x and almost all t, 

{CL) the vector fields f{t, x), gi{t,x), . . . , gm{t, x) are continuous on [0, ti] x X, 

• the vector fields f{t,x), gi{t,x), gm{t,x) have bounded components 
on [0,ti] X X, 

• the function fo{t, x, u) is convex in u for all {t, x) G [0, ti] x X, 

• X^ G Axo{ti). 

Then there exists an optimal control u G L2([0,ti],M"*) for the problem (5.9)- 
(5.12). 
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For Euler's problem embedded into (5.2)-(5.8), we have: 

• m = 1, 

• the set X = $(M) is compact, 

• the function /o = is continuous, nonnegative, and convex. 



d ^ X1X2 d ^ X2 d 
xf + x^dxi xl+xldx2 + x\ dxa' 



the vector fields f{x) = ^ ^2 ^ + ^2 , ^2 ^ + 



, , d d . , , , , , 

gux) = —X2^ \- Xi— — are contmuous and have bounded components 

0x1 0x2 

on X, 

• x^ & Axo{t\) as supposed in (4.3). 

So all hypotheses of Th. 5.1 are satisfied, and there exists optimal control u G 
L2[0,ti] for Euler's problem. 

5.3 Boundedness of optimal controls 

One can prove essential boundedness of optimal control in Euler's elastic prob- 
lem by virtue of the following general result due to A.Sarychev and D.Torres. 

Theorem 5.2 (Th. 1 [35]). Consider an optimal control problem of the form 
(5.9)-(5.12). Let fo S C^i[0,ti] x X x M'",R), /, g, e C^{[0,ti] x X;R"), 
i = 1,. . . ,m, and ip{t,x,u) = f{t,x) + Y,"Li Uigi{t,x). 
Under the hypotheses: 

(HI) full rank condition: dimspan(gi(t, x), . . . ,gm{t,x)) = m for all t G [0,ti] 
and X e X; 

{H2) coercivity: there exists a function 6 : M. ^ M. such that fo{t,x,u) > 
Oi\\u\\)>C y (t,x,u) e[0,ti]x X xR"", and lim 7.^=0; 

{H3) growth condition: there exist constants 7, /3, 77, and ji, with 7 > 0, /3 < 2, 
and n > max{/3 — 2, —2}, such that, for allt G [0, ti], x G X, and u G K™, 
it holds that 

{\fot\+\foxi\ + \\foVt-foM\ + \\foVxi-foxM\)hr <lfo+V, i=l,...,n, 

all optimal controls u(-) of the problem (5.9)-(5.12) which are not abnormal 
extremal controls, are essentially bounded on [0,t-i]. 



It is easy to see that all hypotheses of Th. 5.2 hold: 

9^ + ^13^ 



(HI) g{x) = -X2£- + x^£-^0onX■ 



iH2) 9ir) = r2; 
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{H3) fat = fax, = ^Pl = 0, \\(px,\\ < C on X. The required bound ||/o<Pa;J| • 
ll^ll^ < 7/(f + ?7 is satisfied for fi = 1, /i = 1, 7 = C, ry = 0. 

Thus in Euler's elastic problem all optimal controls which are not abnormal 
extremal controls are essentially bounded: w(-) € L^[0,ti]. In Subsec. 6.2 
we describe abnormal extremal controls, and obtain a similar inclusion for all 
optimal controls. 

Meanwhile we cite one more general result valid for Euler's problem as well. 

Corollary 5.1 (Cor. 1 [35]). Under conditions of Th. 5.2, all minimizers of 
the problem (5.9) -(5.12) satisfy the Pontryagin Maximum Principle. 

We summarize or results for Euler's elastic problem derived in this section. 

Obviously, we can return back from the problem (5.2)-(5.8) in ^^1x2x3 
initial problem (3.4)~(3.9) in Rl^y x S^. 

Theorem 5.3. Let qi € Ag„{ti). 

(1) Then there exists an optimal control for Euler's problem (3.4) -(3.9) in the 

class u{-) e ^2(0, ^i]. 

(2) // the optimal control u(-) is not an abnormal extremal control, then u{-) € 
Loo[0,ti\. The corresponding optimal trajectory g(-) is Lipschitzian. 

(3) All optimal solutions to Euler's problem satisfy the Pontryagin Maximum 
Principle. 

Certainly, Th. 5.3 is not the best possible statement on regularity of solutions 
to Euler's problem. We will derive from Pontryagin Maximum Principle that 
optimal controls and optimal trajectories are analytic, see Th. 6.3. 



6 Extremals 

6.1 Pontryagin Maximum Principle 

In order to apply Pontryagin Maximum Principle (PMP) in invariant form, 
we recall the basic notions of the Hamiltonian formalism [2,21]. Notice that 
the approach and conclusions of this section have much intersection with the 
book [21] by V.Jurdjevic. 

Let M be a smooth n-dimensional manifold, then its cotangent bundle T*M 
is a smooth 2n-dimensional manifold. The canonical projection tt : T*M M 
maps a covector A G T*M to the base point q £ M. The tautological 1-form 
s e A^(T*M) on the cotangent bundle is defined as follows. Let A e T*M and 
V € T\(T*M), then (sx.v) = (A,7r*w) (in coordinates s = pdq). The canonical 
symplectic structure on the cotangent bundle a G A^(T*M) is defined as a = 
ds (in coordinates a = dp A dq). To any Hamiltonian h e C°°{T*M), there 
corresponds a Hamiltonian vector field on the cotangent bundle h € Vec(T*M) 
by the rule a\{-,h) = dxh. 
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Now let M = R^ j^ X S\ be the state space of Euler's problem. Recall that 
the vector fields 

d d d d d 

Xi=cos6— +saie—, ^2=— , X3=siii^- cos6'— 

ox ay oU ox oy 

form a basis in the tangent spaces to M. The Lie brackets of these vector fields 
arc given in (3.10), (3.11). Introduce the linear on fibers of T*M Hamiltonians 
corresponding to these basis vector fields: 

/i,(A) = (A,Xi), XgT*M, i= 1,2,3, 

and the family of Hamiltonian functions 

h^^{X) = (A, Xi + UX2) + = /11(A) + uh2{X) + 

A e T*M, u e M, ugR. 

the control-dependent Hamiltonian of PMP for Euler's problem (3.4)-(3.9). 

By Th. 5.3, all optimal solutions to Euler's problem satisfy Pontryagin Max- 
imum Principle. We write it in the following invariant form. 

Theorem 6.1 (Th. 12.3 [2]). Let u{t) and q{t), t G [0,ti], be an optimal control 
and the corresponding optimal trajectory in Euler's problem (3.4)-(3.9). Then 
there exist a curve At S T*M, n{Xt) = q{t), t e [0,ti], and a number u < for 
which the following conditions hold for almost all t £ [0,ti]: 

At = K{tMt) = hi{Xt) + u{t)h2{Xt), (6.1) 
Kit)i>^t)=ma^h-:,{Xt), (6.2) 

(i^,At)^0. (6.3) 

Using the coordinates {hi, /i2, hs, x, y, 9), we can write the Hamiltonian sys- 
tem of PMP (6.1) as follows: 

hi = -uh3, (6.4) 

/i2 = /13, (6.5) 

/13 = uhi, (6-6) 

x = cos^, (6.7) 

y = sin6', (6.8) 

e = u. (6.9) 

Notice that the subsystem for the vertical coordinates (/ii, /12, /is) (6.4)-(6.6) is 
independent of the horizontal coordinates {x, y, 9), this is a corollary of the left- 
invariant symmetry of system E and of appropriate choice of the coordinates 

{hi,h2,h3), see [2]. 

As usual, the constant parameter v can be either zero (abnormal case), or 
negative (normal case, then one can normalize 1/ = — 1). 
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6.2 Abnormal extremals 



Consider first the abnormal case: let z/ = 0. The maximality condition of 
PMP (6.2) reads: 

hl(X) = hi(X) + u/i2(A) max, (6.10) 

thus h2{Xt) = along an abnormal extremal At. Then Eq. (6.5) yields /i3(At) = 
0, and Eq. (6.6) gives u(t)hi{Xt) = 0. But in view of the nontriviality condition 
of PMP (6.3), we have hi{Xt) ^ 0, thus u{t) = 0. So abnormal extremal 
controls in Euler's problem are identically zero. Notice that these controls are 
singular since they are not uniquely determined by the maximality condition of 
PMP (6.10). 

Now we find the abnormal extremal trajectories. For u = the horizontal 
equations (6.7)-(6.9) read 

{X = cos 9, 
y = sme, 
6 = 0, 

and the initial condition {x,y,6){0) = (0,0,0) gives 

x{t) = t, y{t) = 0, e{t) = 0. 

The abnormal extremal trajectory through qo = Id is the one-parameter sub- 
group of the Lie group E(2) corresponding to the left-invariant field Xi. It is 
projected to the straight line {x,y) = {t,0) in the plane {x,y). The correspond- 
ing elastica is a straight line segment — the elastic rod without any external 
forces applied. This is the trajectory connecting go to the only attainable point 
qi at the boundary of the attainable set Aqg{ti). 

For u = the elastic energy is J = 0, the absolute minimum. So the 
abnormal extremal trajectory q(t), t € [0,ti], is optimal; it gives an optimal 
solution for the boundary conditions qo = (0,0,0), qi = (ti,0,0). 

Combining the description of abnormal controls just obtained with Th. 5.3, 
we get the following statement. 

Theorem 6.2. For any qi G Aqg{ti), the corresponding optimal control for 
Euler's problem (3.4) -(3.9) is essentially bounded. 

6.3 Normal case 

Now let u = —1. The maximality condition of PMP (6.2) reads 

= hi + uh2 — liu"^ max, 
2 «eM 

whence = h2 — u = and 

ou 

u = h2. (6.11) 
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The corresponding normal Hamiltonian of PMP is 

and the normal Hamiltonian system of PMP reads 



A = H{\) 



h2 = hs, 
/13 = hihs, 
^q = Xi + /i2^2- 



(6.12) 



This system is analytic, so we obtain the following statement (taking into ac- 
count analyticity in the abnormal case, see Subsec. 6.2). 

Theorem 6.3. All extremal (in particular, optimal) controls and trajectories 
in Euler's problem are real-analytic. 

Notice that the vertical subsystem of the Hamiltonian system (6.12) admits a 
particular solution {hi, /12, ^3) = (0,0,0) with the corresponding normal control 
w = /i2 = 0. Thus abnormal extremal trajectories are simultaneously normal, 
i.e., they are not strictly abnormal. 

We define the exponential mapping for the problem: 

Exp,^ : T;M ^ M, Exp,^(Ao) = tt o e'^"{Xo) = q{ti). 
The vertical subsystem of system (6.12) has an obvious integral: 



hi + h 



const > 0, 



and it is natural to introduce the polar coordinates 

hi = r cos a, hz = r sin a. 
Then the normal Hamiltonian system (6.12) takes the following form: 

a = /12, 
/i.2 = rsina, 
r = 0, 
X = cos 6, 
y = sin 9, 
j = h2. 

The vertical subsystem of the Hamiltonian system (6.13) reduces to the 
equation 

d = r sin a. 



(6.13) 



In the coordinates 



/l2, 



/? = a + TT 



23 



we obtain the equation of pendulum 

f/3 = c, 

1 c = — r sin /? 

known as Kirchhoff's kinetic analogue of elasticae. Notice the physical meaning 
of the constant: 

r=|, (6.14) 

where g is the gravitational acceleration and L is the length of the suspension 
of the pendulum, see Fig. 6. 




Figure 6: 
Pendulum 



The total energy of the pendulum is 



E. 



■ rcos/3 e [— r, +cxd), 



(6.15) 



this is just the Hamiltonian H. 

The geometry of solutions to the vertical subsystem of the normal Hamilto- 
nian system (6.12) 



hi = /i3, 
/i3 = hihs 



I3 = c, 

c = — r sin/J, 
f = 



is visualized by intersections of level surfaces of the integrals 



(6.16) 



H=hi 



hi + h 



Depending upon the structure of intersection of circular cylinders = const 
with parabolic cylinders H = const, the following cases are possible, see Figs. 7- 
12: 



(a) H = -r,r>0 
(0,0), 



pendulum stays at the stable equilibrium (/?, c) 
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(b) H E {—r,r), r > ^ pendulum oscillates between extremal values 
of its angle f3, 

(c) H — r, r > pendulum either stays at the unstable equilibrium 
(/?, c) = (tt, 0) or tends to it for infinite time, 

(d) H > r, r > pendulum rotates non- uniformly counterclockwise 
(ft.2 > 0) or clockwise (/i2 < 0), 

(e) H = r = => pendulum is immovable in the absence of gravity, with 
zero angular velocity h2, 

(J) H > r ~ =^ pendulum rotates uniformly in the absence of gravity 
/i2 counterclockwise (/12 > 0) or clockwise (/12 < 0). 



Figure 7: H = -r , r > Figure 8: H £ {-r,r), r > 

The equation of mathematical pendulum is integrable in elliptic functions. 
Consequently, the whole Hamiltonian system (6.13) is integrable in quadratures 
(one can integrate first the vertical subsystem, then the equation for 9, and then 
the equations for x, y). In Sec. 8 we find explicit parametrization of the normal 
extremals by Jacobi's elliptic functions in terms of natural coordinates in the 
phase space of pendulum (6.16). 

First we apply continuous symmetries of the problem. The normal Hamil- 
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Figure 9: H = r > 




tonian vector field reads 

d d 
H = -h2h3—j- + h3—j- 
ohi a/i2 



dhs 



dx 



dy 



■89 



, 9 . d ^ d . ^ d , d 

= h rsina— — + cos0— +sin6'— + /i2^- 

oa o/i2 ox ay 06 



The Hamiltonian system (6.13) is preserved by the one-parameter group of 
transformations 



(a, r, h2,x, y, 6, t) {a, re h2e xe^, ye*, 9, te^) 



(6.17) 



obtained by continuation to the vertical coordinates of the group of dilations of 
the plane Rl y (3.13). 

The one-paxameter group (6.17) is generated by the vector field 

„ ^ d , d d d 

Z = -^r— -h2jrj-+ X— + y — . 
or on2 ox ay 

We have the Lie bracket and Lie derivatives 



[Z,H] = -H, 

Zr = -2r, Zh2 = -^2, 



Hr = 0, Hh2 = rsina. 



(6.18) 
(6.19) 



The infinitesimal symmetry Z of the Hamiltonian field H integrates to the 
symmetry at the level of fiows: 



„t'H „ sZf 



oe«^(A) = e*^oe*-^(A), t' = eH, 



A e T*M. 



The following decomposition of the preimage of the exponential mapping N into 
invariant subsets of the fields H and Z will be very important in the sequel: 





7 

= iv = y iVi, 

i=l 




(6.20) 


iVi = 


{\eN\r^Q, Ee{- 


r)}, 


(6.21) 


N2 = 


{Xe N \ ry^O, E€{r 


,+^)} = N+UN^, 


(6.22) 


Ns = 


{XeN \ ry^O, E = r, 




(6.23) 


iV4 = 


{XgN \ ry^O, E = - 


r}, 


(6.24) 


A^5 = 


{X&N\r^Q, E = r, 




(6.25) 


Ne = 


{Ae Ar|r = 0, c^O} 


= N+UN^, 


(6.26) 


N-j = 


{A e AT 1 r = c = 0}, 




(6.27) 


Nt = 


-- NiU{Xe N \ sgnc = 


±1}, i = 2, 3, 6. 


(6.28) 



Any cylinder {X G N \ r = const ^0} can be transformed to the cylinder 
C = {X & N \ r = 1} hy dilation Z; the corresponding decomposition of the 
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phase space of the standard pendulum 

\c=—smp, 
is shown at Fig. 14, where 

a = N,n{r=l}, z = 1, . . . , 5. 





Figure 13: Phase portrait of pendu- 
lum 



Figure 14: Decomposition of the 
phase cylinder of pendulum 



In order to integrate the normal Hamiltonian system 

'/3 = c, 
c = — rsin/3, 

X = H{X) ^ K = 0' (g_29) 

X — cos 0, 

y = sin 9, 

e = c, 

we consider natural coordinates in the phase space of the pendulum. 

7 Elliptic coordinates 

The twenty-sixth of the "Lectures on dynamics" by C.Jacobi [16] is called "El- 
liptic coordinates" and begins by the well-known words: "The main problem 
in the integration of these differential equations is the introduction of conve- 
nient variables, there being no general rule for finding them. Therefore one 
has to adopt the opposite approach and, finding a remarkable substitution, to 
seek the problems for which this substitution can be successfully used". Note 
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that the coordinates introduced below are unrelated to Jacobi's elliptic coor- 
dinates. Moreover, our procedure was opposite to that described by Jacobi: 
we introduced our elliptic coordinates specifically of parameterizing extremals 
and finding Maxwell points in generalized Dido's problem [30-33] and in Euler's 
problem. Elliptic coordinates lift the veil of complexity over the problems gov- 
erned by the pendulum equation and open their solution to our eyes (see Fig. 15). 
Here we have an important intersection point with Jacobi: our coordinates are 
introduced by using Jacobi's elliptic functions, see Sec. 13 and [22], [38]. An- 
other important moment will be the study of conjugate points, that is, solutions 
to Jacobi equation, along extremals [34]. 



c 




Figure 15: Elliptic coordinates in the phase 
cylinder of pendulum 



7.1 Time of motion of the pendulum 

Elliptic coordinates in the phase cylinder of the standard pendulum 

l'^^"'. ^ (Ac) eC^^i xM„ (7.1) 

c = — sm p. 
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were introduced in [30] for integration and study of the nilpotent sub-Riemannian 
problem with the growth vector (2,3,5). Here we propose a more natural and 
efficient construction of these coordinates. 
Denote 

P = M+cXMt, C = CiUC^UC^, 
and consider the mapping 

$ : P^C, $ : (c,i)^(A,Ct), 

where {(3t,Ct) is the solution to the equation of pendulum (7.1) with the initial 
condition 

Po = 0, Co = c. (7.2) 

The mapping $ : P ^ is real-analytic since the equation of pendu- 
lum (7.1) is a real-analytic ODE. 

First we show that $ is a local diffeomorphism, i.e., the Jacobian 



d{Pt,ct) 



d{c,t) 

By virtue of system (7.1), we have 

Ct, 



df3t dl3t 

act oct 
dc dt 



^ V(c, t) e P. 



cm 

dt 



dct 

dt 



Further, denote 



Since 



dc 



dct 
dc 



■sin/?t. 







d dPt 


dc^ 


dt 


Kdcj 


dc dt 


dc 






_ d dct _ 


dc^ 


dt 


m 


~ dc dt ~ 



vit), 



■cos/3tC(i), 



the pair (((t),r]{t)) is the solution to the Cauchy problem 

dPo 



C = v, 

?7 = - cos (3t C, 

Now consider the determinant 

diPuCt) 



C(0) = 
77(0) = 



dc 
dcp 
dc 



0, 
1. 



d{t) 



d{c,t) 



m 



Ct 

■sin/3t 



(7.3) 
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Differentiating by rows, we obtain 

d{t) 



vit) 



■ sin Pt 
■sinPt 



thus 



d{t) = d{0) 



c 

1 



at) 

■cosACW -cosACW 



(7.4) 



Denote by 
k2 



E + 1 
2 







sm 



sin^ f + f 



e (0,1), 



(Ac) GCi, 

(/3,c) eC2, 



a reparametrized energy (iJ = — cos/3) of the standard pendulum; below ki, 
k2 will play the role of the modulus for Jacobi's elliptic functions, and 

r-^^ dt 



K{k) = / 
Jo 



Vl - fc2 sin^t' 



fee (0,1), 



the complete elliptic integral of the first kind, see Sec. 13 and [22]. It is well 
known [22] that the standard pendulum (7.1) has the following period of motion 
T depending on its energy E: 

-KE <1 4^ (/3, c) e Ci ^ T = 4K{ki), (7.5) 
E=l, Py^TT ^ (/3,c)eC3 => T = oo, (7.6) 
E>1 ^ (/3,c)eC2 => T = 2K{k2)k2. (7.7) 

Introduce the equivalence relation ~ in the domain P as follows. For (ci, fi) € P, 
(c2,i2) € P, we set (ci,<i) ~ (02,^2) iff ci = C2 = c and 



i2 = ti (mod 4if(fci)), fci = - 

2 

i2 = ti (mod 2K{k2)k2), ^2 = - 



for ce (0,2), 

for c = 2, 

for c e (2, +00). 



That is, we identify the points {c,t\), {c,t2) iff the corresponding solutions to 

the equation of pcnduhim (7.1) with the initial condition (7.2) give the same 
point {Pt, Ct) in the phase cylinder of the pendulum 5*^ x at the instants ti, 

Denote the quotient P = P/ ^. In view of the periodicity properties (7.5)- 
(7.7) of the pendulum (7.1), the mapping 



$ : P ^ C, 



$(c,t) = {f3t,ct), 
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is a global analytic difFeomorphism. Thus there exists the inverse mapping, an 
analytic diffeomorphism 



F : C ^ P, 
Fi^,c) = ico,ip), 



(7.8) 



where tp is the time of motion of the pendulum in the reverse time from a point 
(/?, c) G C to the semi-axis {/3 = 0, c = 0}. In the domains Ci and C^, the 
time is defined modulo the period of the pendulum 4:K{ki) and 2K{k2)k2 
respectively. 

We summarize the above construction in the following proposition. 

Theorem 7.1. There is an analytic multi-valued function 

<fi : d = CiUC}uC^ 

such that for Pq = 0, cq > and the corresponding solution {l3t,Ct) of the Cauchy 
problem (7.1), (7.2), there holds the equality 



t (mod AK{ki)) 
t 

t (mod 2K{k2)k2) 



for (A,cO eCi, 
for (A,q) €^2+, 
for (/3„c,) eC+. 



In other words, (^(/3f , ct) is the time of motion of the pendulum in the reverse 
time from the point {(3t,ct) € C to the semi-axis {/? = 0, c > 0}. 

7.2 Elliptic coordinates in the phase space of pendulum 

In the domain Ci U C2 U C3, wc introduce elliptic coordinates {(p, k), where ip is 
the time of motion of the pendulum from the semi-axis {/3 = 0, c > 0} (in the 
domain C = Ci UC^ UC^) or from the semi-axis {/? = 0, c < 0} (in the domain 
C = C2 U C^), and k G (0, 1) is a reparametrized energy of the pendulum — 
the modulus of Jacobi's elliptic functions. 



7.2.1 Elliptic coordinates in Ci 

If (/?, c) e Ci, then we set 

sin ^ = ki sniip. ki), 

- = fci Cn((y9,fcl), 

cos ^ = dn(<^, fci). 



ki 



ip (mod 4K{ki)) € [0,4K{ki)]. 



(7.9) 



(7.10) 
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Here and below cn, sn, dn are Jacobi's elliptic functions, see Sec. 13 and [22]. 

The function ^p thus defined is indeed the time of motion of the pendulum 
from the semi-axis {/3 = 0, c>0}in view of the following: 

(/3 = 0, OO) ^ ip = Q, (7.11) 



dt 



= 1, (7.12) 

Eq. (7.1) 



the total derivative w.r.t. the equation of pendulum (7.1). 

The mapping {13. c) (fc, ip) is an analytic diffeomorphism since it decom- 
poses into the chain of analytic diffeomorphisms: 

(/3,c)^i'(co,v')^^(fci,^), 
where (a) is defined by F (7.8), while (6) is given by 

l E + l Co 

compare with (7.10). 



7.2.2 Elliptic coordinates in 

Let (/3, c) G C2 ■ Elliptic coordinates {ip,ki) in the domain C2 are analytic 
functions (<p, fci) defined as follows: ip is the time of motion of the pendulum 
from the semi-axis {/? = 0, c > 0}, and fci = By the uniqueness theorem 

for analytic functions, in the domain C2 we have the same formulas as in Ci: 

sin^ = fcisn(.^, fci), (7.13) 
I = fcicn(^,A;i), (7.14) 

cos ^ = dn(y, fci), (7.15) 
fci = y^±^G (l,+oo). 

Here Jacobi's elliptic functions sn(u, fci), cn(u, fci), dn(w, fci) for the modulus 
ki> \ are obtained from those defined in (13.1)-(13.3) by the analytic continu- 
ation along the complex modulus /ci e C through the complex plane around the 
singularity fci = 1, see Sec. 3.9 and Sec. 8.14 [22]. In order to obtain Jacobi's 
functions with the modulus in the interval (0, 1), we apply the transformation 
of modulus ^, see formulas (13.7), (13.8) in Sec. 13. Transforming equali- 
ties (7.13)-(7.15) via formulas (13.7), (13.8), we obtain the following expressions 
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for elliptic coordinates (</?, ^2): 




(7.16) 



Certainly, one can verify directly that <^ is indeed the time of motion of the 
standard pendulum from the point c) to the semi-axis {/3 = 0, = c > 0} in 
the reverse time by checking the conditions (7.11), (7.12) in the domain C^, 
but our idea is to obtain "for free" equalities in Cj" from equalities in Ci via 
the transformation of the modulus fc 1— > p 



7.2.3 Elliptic coordinates in 

Let (/?, c) G C^. Elliptic coordinated on the set are given by {ip,k = 
1), where if is the time of motion of the pendulum from the semi-axis {/3 = 

0, c > 0}; and k = \J^^ = 1- The analytic; expressions for Lp are obtained 
by passing to the limit ki 1 — in formulas (7.9) or to the limit k2 — > 
1 — in formulas (7.16), with the use of formulas of degeneration of elliptic 
functions (13.6). As a result of the both limit passages, we obtain the following 
expression for the elliptic coordinate (p on the set C^: 



sin ^ = tanh ip, 

c _ 1 

2 cosh tfi ' 

P 1 
cos — = — , — . 
2 cosh 



7.2.4 Elliptic coordinates in U C3 

For a point (/3, c) G C = U , elliptic coordinates {cp, k) cannot be defined 
in the same way as in C = Ci U C2 U since such a point is not attainable 
along the flow of the pendulum (7.1) from the semi-axis {/3 = 0, c > 0}, see the 
phase portrait at Fig. 13. Now we take the initial semi-axis {/3 = 0, c < 0}, and 
define inC equal to the time of motion of the pendulum from this semi-axis 
to the current point. That is, for points (/?, c) € (7 we consider the mapping 

F(c,t) = (A,Ct), c<-2, 
/3o = 0, Co = c. 
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and construct the inverse mapping 

$(/3,c) = (co,<p). 

The pendulum (7.1) has an obvious symmetry — reflection in the origin 
(/3 = 0, c = 0): 

i : (/3,c)^(-/?,-c). (7.17) 



In view of this symmetry, we obtain: 

$(/3,c) = (co,</?), 

$(-/?, -C) = {-Co, if), 



(/3,c)GC2-UCf, 
(-/3,-c) eC2+UC3+, 



thus 



ifiiP, c) = ^(-/3, -c), (/?, c) G U C3- 



On the other hand, the energy of the pendulum E and the modulus of elliptic 
functions ^2 are preserved by the reflection (7.17). So we have the following 
formulas for elliptic functions in C. 



sm- = - 
c 1 



- sn 



dn 



(/3,c)eC3- 



sm 

c 
2 

cos 



f3 
2 



^2' 



cos ^ = cn ( A:2 ) , 



- tanh (fi, 
1 



cosh (j3 ' 



2 cosh (f 



Summing up, in the domain Ci U C2 U C3 the elliptic coordinates (</?, A:) are 
deflned as follows: 



(Ac) eCi 



£; + 1 



sin^ = kisn{ip,ki), 
- = kicn{ip,ki), 
cos ^ = dn{ip,ki), 



e (0,1), (mod AK{ki)) G [0,4i4:(fci)], 
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= y^g— ^ e (0,1), if (mod 2i^(fc2)fe) e [0,2if(fc2)A;2], ± = sgnc, 



sin — — ± tanh <p, 
2 

- = ± — , 
2 cosh (/? ' 

/3 1 
cos - = — —, 
2 cosh if 

± = sgn c. 



Remark. In such a definition of elUptic coordinates, the domains C2 and 
(C^ and C^) have different status: the coordinate ip is discontinuous when 
crossing the scparatrix , and analytic on the separatrix . This is a c;onse- 
quence of the fact that in Ci, C2 , the coordinate tp is defined uniformly — 
as the time of motion of the pendulum from the semi-axis {/3 = 0, c > 0}, while 
in C2 , this is the time of motion from another semi-axis {/3 = 0, c < 0}. 
This different status is refiected in the fact that elliptic coordinates in C2 are 
obtained from elliptic coordinated in C\ by analytic continuation (with the use 
of the transformation k 1-^ j: of Jacobi's functions), after which elliptic coordi- 
nates in C2 are obtained from C2 via the symmetry i of the pendulum (7.17). 

The use of analytic continuation from Ci to Cj" allows us to obtain "gratis" 
all formulas in C2 from the correspoiidiiig formulas in Ci via the transforma- 
tion k (13.7), (13.8). As usual for analytic functions, analytic continuation 
respects only equaUties; inequalities are not continued in such a way, in partic- 
ular, we will have to obtain bounds for roots of equations independently in Ci 
and C2. But in order to obtain equalities in the cylinder C (and in the preimage 
of the exponential mapping N) we will make use of the following chain: 



C4 Ci — 



(7.18) 



fe^l-O 


i 


1 ^ 


Y 



Such a chain will be useful not only in Euler's problem, but in all problems 
governed by the pendulum (7.1), e.g. in the nilpotent (2,3,5) sub-Riemannian 
problem [30-33], in the plate-ball problem [21], in the sub-Riemannian problem 
on the group of motions of the plane. 
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At Fig. 15 wc present the grid of elliptic coordinates in the phase cylinder of 
the standard pendulum (M.^ ^ x 5^). In the domain C\ (oscillations of pendulum 
with low energy E <1) we plot the curves k = const, = const; in the domain 
C2 (rotations of pendulum with high energy E > 1) we plot the curves k = const, 
ip = const; these domains are separated by the set C2 (motions of pendulum 
with critical energy E = 1), consisting of two separatrices k = 1 and the unstable 
equilibrium. 

7.3 Elliptic coordinates 

in the preimage of the exponential mapping 

In the domain TV = A^i U U JV3 (recall decomposition (6.20) (6.28)), the verti- 
cal subsystem of the Hamiltonian system (6.12) has the form of the generalized 
pendulum 

'P = c, 

< c=-r sin (3, (7.19) 
f = 0. 

V 

Elliptic coordinates in the domain N have the form {ip,k,r). On the set Ni U 
A^2^ U , the coordinate ip is equal to the time of motion of the generalized 
pendulum (7.19) from a point (/3 = 0, c = cq > 0, r) to a point (/3, c, r), while on 
the set U the time of motion is taken from a point (/3 = 0, c = Cq < 0, r). 
The one-parameter group of symmetries 

{0,c,r,t)^{P,ce-\re-^\te') 

of the generalized pendulum (7.19) is a restriction of action of the group (6.17). 

We apply this group to transform the generalized pendulum (7.19) in the domain 
{r > 0} to the standard pendulum (7.1) for r = 1. This transformation preserves 
the integral of the generalized pendulum 

lE + r 1.2/3 
V 2r V 2 4r 

Thus we obtain the following expressions for elliptic coordinates in the do- 
main N from similar expressions in the domain C, see Subsec. 7.2. 

sin^ = ki sn{y/rip,ki), 
< 2 = ki^/rcii{^/r(p, ki), 
cos ^ = dn(v^</j, fci), 

y/np (mod AK{ki)) e [0,AK{ki)], 



X={(3,c,r)em ^ 
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X = {f3,c,r)eNi 



k2 



2r 
E + r 



G (0,1), 



Vftp (mod 2K{k2)k2) G [0,2i^(fc2)/c2] 



± = sgn c. 



A=(/3,c,r)e7V3± 



1, 



V G 



sin — = ±tanh(\/r(/j), 



2 cosh(Y^</?) ■ 
/3 1 



cos 



2 cosh{y/rip) ' 



± = sgn c. 



In the domain it will also be convenient to use the coordinates {k2, ip, r), 
where 

V> = ^, (mod 2i<:(A;2)) G [0, 2ii:(A;2)]. 

fc2 

In computations, if this does not lead to ambiguity, we denote the both moduli 
of Jacobi's functions ki and k2 by k, notice that k G (0, 1), this is the normal 
case in the theory of Jacobi's elliptic functions, see [22]. 



8 Integration of the normal Hamiltonian system 
8.1 Integration of the vertical subsystem 

In the elliptic coordinates {(f, k, r) in the domain N , the vertical subsystem (7.19) 
of the normal Hamiltonian system A = H{\) rectifies: 

</9 = 1, = 0, r = 0, 

thus it has solutions 

ipt = if + t, k = const, r = const . 

Then expressions for the vertical coordinates (/?, c, r) are immediately given by 
the formulas for elliptic coordinates derived in Subsec. 7.3. For X G N \ N, 
the vertical subsystem degenerates and is easily integrated. So we obtain the 
following description of the solution (/?t,Ct,r) to the vertical subsystem (7.19) 
with the initial condition {Pt,Ct,r)\^^Q = {l3,c,r). 
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A e 



A e 



XGN4 
AGiVs 
A e iVe 



sin ^ = ki sn{\/ript), 
cos y = dn(\/''¥'t), 

sm — = ± sn — ; — 
2 \ k 

cos-=cn(^ — 
2 k \ k 



( . Pt 

sm 



cos 



±tanh(\/r<^t), 
1 

2 coslifN/r^Sf ) ' 

^ ^ _^ V'l' 

2 



"cosh(i/r(pt) ■ 
A = 0, ct = 0. 

/3t = TT, Ct = 0. 

A = + /3, Ct = c. 
Ci = 0, r = 0. 



8.2 Integration of the horizontal subsystem 

The Cauchy problem for the horizontal variables (a;, y, 6) of the normal Hamil- 
tonian system (6.29) has the form 



X = cos = 2 cos^ 



e 



w = sin ^ = 2 sin - cos - , 
" 2 2' 



Xo 


= 0, 


yo 


= 0, 


Oo 


= 0, 



thus 



= /3t - /?. 



(8.1) 



We apply known formulas for integrals of Jacobi's elliptic functions, see Sec. 13, 
and obtain the following parametrization of normal extremal trajectories. 
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If Ae7Vi,then 

Q 

sin — = kd\i{^/rLp)sn{^/ript) — ksn-{^/rLp) &si{\frLpt), 

Q 

COS — = dn{y/rip) dn{y/ript) + k"^ sn{y/r(f) sn{y/ript), 
xt = ^ dn'(^/f^)(E(^/f<^^) - E(^/^^)) 

+ dn(v^(/?) sn(v^(^)(cn v^iyj) — cn(v^(pt)) 
2)^2 

+ -j= sn\^^)i^t + Ei^ip) - Ei^ift)) - t, 
2k 

yt = ^(2dn^(v^<p) - l)(cn(Vr-(yj) - cn{./ript)) 
^/r 

2k 

- — An{^^){2{^{^^t) - E(x/F(^)) - ^prt). 

Here E(u, A:) is Jacobi's epsilon function, see Sec. 13 and [22]. 

The parametrization of trajectories in A^2^ is obtained from the above parametriza- 
tion in A^i via the transformation fc i— s- j; described in Subsubsec. 7.2.4; af- 
ter that, trajectories in N2 arc obtained via the reflection i (7.17), see the 
chain (7.18). In the domain N2, we will use the coordinate 



Then we obtain the following. 
If A e Nf, then 



sin — = ±{ci\{^ij) ani^i/jt) - sni^/rtp) cn{y/ript)) , 
cos = cn{y/rip) cn{y/ript) + sn^y/rij}) sn{y/ript), 

Xt = ^(1 - 2sn2(^V)) (^(E(^^^V^) - E(^A^V)) - ^^^^*) 
4 

+ - — ^ cn(^ip) SIl(^'^jJ)(dIl(^/rtp) — dn( ^/r'^jJt)), 
k^/r 

yt = ± (^^(2cn2(^V) - l)(dn(^/rV) - dn(^/r^t)) 

2 / 2 2- k^ 

--^sn(\/^V')cn(v^V) U(E(v^Vt) -E(v^^)) - —j^V^t 

The formulas in are obtained from the above formulas in N2 via the 
limit A; ^ 1 — 0, see the formulas of degeneration of Jacobi's functions (13.6), 
and compare with chain (7.18). 
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Consequently, if A e N^, then 

^ _ _l_ / tanh(v^yt) _ tanhA/ry)) \ 

2 V cosh{^/r^p) cosh{^/r(fit) J ' 

9 1 
cos ^ = — + tanh{^/r(p) ta.nh{^/npt), 

2 1 rz. \\M , 4tanh(-;/r(^) 



a;j = (1 - 2tanh'^(yf(p))t + 



yt = ± ^ —2-^. 1 



■v/rcosh(y^<^) \cosh(y^(^) COsh(y^(^t] 

1 1 



\/r Vcosh^ ^(^) y Vcosh(^(^) cosh(v^(pt) 
g tanh(Vry) ^ 



COSh(y^yi) 

Now we consider the special cases. 
If A e iV4 U iV5 U A^T, then 

Qt = 0, a;t = Vt = 0. 

If A e A/e, then 

sin ct 1 — cos ct 



= ct, xt = , Vt = 



c c 
So we parametrized the exponential mapping of Euler's elastic problem 

Expj : X={p,c,r)^qt = {et,xt,yt), XeN = T*^M, qt&M, 

by Jacobi's elliptic functions. 



8.3 Euler elasticae 

Projections of extremal trajectories to the plane (x, y) are stationary configura- 
tions of the elastic rod in the plane — Euler elasticae. These curves satisfy the 
system of ODEs 

X = cos 0, 
y = sinO, 

e = -rsm{e-p), r,/3 = const. (8.2) 

Depending on the value of energy ^ ~ ~^ — rcos(0 — /?) G [— r, +oo) and the 

constants of motion r G [0, +00), (3 G S^, of the generalized pendulum (8.2), 
elasticae have different forms discovered by Euler. 

If the energy E takes the absolute minimum —r 7^ 0, thus A G N4, then the 
corresponding elastica {xt,yt) is a straight line (Fig. 16). The corresponding 
motion of the generalized pendulum (Kirchoff 's kinetic analogue) is the stable 
equilibrium. 
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Figure 16: E = ±r, r > 0, c = Figure 17: E e (-r,r), r > 0, G 




Figure 18: E e (-r,r), r > 0, k = 
1 




Figure 20: e (-r,r), r > 0, fc = 




Figure 19: E e (-r, r), r > 0, fc e 




Figure 21: E G (-r,r), r > 0, fc e 
(fco.l) 
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Figure 22: E = r>0, /3^tt 



Figure 23: E>r>0 




li E E (— r, r), r 7^ 0, thus \ E Ni, then the pendulum oscillates between 
extremal values of the angle, and the angular velocity 9 changes its sign. The 
corresponding elasticae have inflections at the points where ^ = 0, and vertices 
at the points where \9\ = max since 6 is the curvature of an elastica (.x^, jji). Such 
elasticae are called inflectional. See the plots of different classes of inflectional 
elasticae at Figs. 17-21. The correspondence between the values of the modulus 

l~E~-\-~r 

of elliptic functions k = \ — - — G (0,1) and these figures is as follows: 

V 2r 

k e (0, Fig. 17, 



V2 

k=^ ^ Fig. 18, 

ke(^,ko) Fig. 19, 



, /2 

k=ko Fig. 20, 

k e (fco, 1) => Fig. 21. 

The value k = ^ corresponds to the rectangular elastica studied by James 
Bernoulli (see Sec. 2). The value ko w 0.909 corresponds to the periodic elas- 
tica in the form of figure 8 and is described below in Propos. 11.5. As it was 
mentioned by Euler, when A; ^ 0, the infiectional elasticae tend to sinusoids. 
The corresponding Kirchhoff's kinetic analogue is provided by the harmonic 
oscillator 9 = —r{9 — /3). 

U E = r and 9 — P ^ tt, thus X E N3, then the pendulum approaches 
its unstable equilibrium [9 — P = ^,6 = 0) along the saddle separatrix, and the 
corresponding critical elastica has one loop, see Fig. 22. 

If i? = r ^ and 9 — P = n, thus A G A^5, then the pendulum stays at 
its unstable equilibrium (6 — /3 = ^,6 = 0), and the elastica is a straight line 
(Fig. 16). 

li E > r 0, thus A G N2, then the Kirchhoff's kinetic analogue is the 
pendulum rotating counterclockwise (^ > X G N2) or clockwise {6 < 

<^ A G iV2~). The corresponding elasticae have nonvanishing curvature 9, 
thus they have no infiection points and are called non-infiectional, see Fig. 23. 
The points where |^| has local maxima or minima are vertices of infiectional 
elasticae. 

If r = and 9^0, thus A G Nq, then the pendulum rotates uniformly: 
one may think that the gravitational acceleration is g = (see the physical 
meaning of the constant r (6.14)), while the angular velocity 6 is nonzero. The 
corresponding elastica is a circle, see Fig. 24. 

Finally, if r = and ^ = 0, thus A G then the pendulum is stationary 
(no gravity with zero angular velocity 9), and the elastica is a straight line, see 
Fig. 16. 

Notice that the plots of elasticae at Figs. 17-23 do not preserve the real ratio 
^ for the sake of saving space. 
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9 Discrete symmetries of Euler's problem 

In this section we lift discrete symmetries of the standard pendulum (7.1) to 
discrete pendulum of the normal Hamiltonian system 



'/?= 


c, 


c = 


— r sin 


f = 


0, 


e = 


c, 


X = 


cos^, 


.y = 


sin^. 



9.1 Reflections in the phase cylinder 
of standeird pendulum 

It is obvious that the following reflections of the phase cylinder of the standard 
pendulum C = 5*^ x Mc preserve the field of directions (although, not the vector 
field) determined by the ODE of the standard pendulum (7.1): 

: (/3,c)^(/3,-c), 

£2 : (/3,c)^(-/3,c), 

£3 : (/3,c)^(-/3,-c), 

see Fig. 25. 




Figure 25: Refiections in the phase Figure 26: Refiections of trajecto- 

cylinder of pendulum ries of pendulum 

These refiections generate the dihedral group — the group of symmetries of 
the rectangle 

i?2 = {Id,£^£^£3} 
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with the multipUcation table 











el 


Id 










Id 










Id 



Notice that the reflections e^, reverse direction of time on trajectories 
of the pendulum, while e"^ preserves the direction of time (in fact, is the 
inversion i defined in (7.17)). 

All reflections preserve the energy of the pendulum E = — — cos 0. 

9.2 Reflections of trajectories of the standard pendulum 

We can define the action of refiections on trajectories of the standard pendulum 
as follows: 

: 7. = m,c,) I s G [0,t]} ^Jl = {(/?:, 4) I s G [0,t]}, 

where 



(/32,c2) = (-A-.,Ct-a), 

(/33, a = (-/?„ -c). 



(9.2) 
(9.3) 
(9.4) 



see Fig. 26. 

All reflections e' map trajectories 7^ to trajectories 7*; they preserve both 



the total time of motion t and the energy E ■ 



cos p. 



9.3 Reflections of trajectories 

of the generaUzed pendulum 

The action of reflections is obviously continued to trajectories of the gener- 
alized pendulum (7.19) — the vertical subsystem of the normal Hamiltonian 
system (9.1) as follows: 



^)\sG[0,t]} 



)\se[o,t]}, 



1,2,3, 



(9.5) 



where the functions c* are given by (9.2) (9.4). Then the reflections 
preserve both the total time of motion t, the energy of the generalized pendulum 

E = — — r cos /?, and the elastic energy of the rod 



1 /■* 

2 Jo 



' 2 



c„ ds. 



46 



(9.6) 



9.4 Reflections of normal extremals 

Now we define action of the reflections on the normal extremals 

Xs = e'"iXo)GT*M, sG[0,t], 
i.e., solutions to the Hamiltonian system 

Cs = -r sin /3s, 
r = 0, 

Js = Xi{qs) + CsX2{qs), 
as follows: 

e' : {A, I s e [0,t]} ^ {Xl \ s e [0,t]}, i = 1,2,3, (9.7) 
Xs = {v,,q,) = {Ps,Cs.r,qs), Xi = {ulqi) = {f3l,c^,r,ql). (9.8) 

Here A* is a solution to the Hamiltonian system (9.6), and the action on the 
vertical coordinates 

: K = (/3„ c„ r)}.-K = (/?:, 4, r)} 

was defined in Subsec. 9.3. The action of reflections on horizontal coordinates 
{0, X, y) is described in the following subsection. 

9.5 Reflections of Euler elasticae 

Here we describe the action of reflections on the normal extremal trajectories 

: {q, = {e,,x,,y,) \ s e [0,t]} ^ {ql = {ei,xi,yi) \ s e [0,t]}. 

Proposition 9.1. Let Xg = {ps,Cs,r,qs) and XI = e*(As) = i(3l,cl,r,ql), s e 
[0,t], be normal extremals defined in (9.7), (9.8). Then the following equalities 
hold: 

{l)6l = 6,..-e„ f1) = f^°^^* sin^, Wx,-x,_. \ 
KVs J \ -smOt cosOt J \ yt-yt-s J 

cos 0t — sin 9t \ ( Xt — Xt-. 



(2) 91 -Ot- Ot.s, [ yl J- \ sine, cose, J \ - 



(3) 



3 

\ I Xg 



Proof. We prove only the formulas in item (1), the next two items are studied 
similarly. By virtue of (8.1) and (9.2), wc have: 

— f^l ~ f^o — l^t-s — Pt = 9t-s — 9t- 
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Further, 



and similarly 



xl = cos 9].dr = I cos{9t-r — Ot) dr 
Jo Jo 

= cos 9t / cos 6t-r dr + sin 9t / sin 9t-r dr 
Jo Jo 

= cos6t{xt - xt-s) + sm.6t{yt - Vt-s), 

yl = sin dr = sin{9t-r — 9t) dr 
Jo Jo 

= cos9t{yt - yt-s) - sm9t{xt - Xts)- 



□ 

Remark. Notice the visual meaning of the action of the reflections on elastica 
{{xs,ys) I s e [0,t]} in the case {xt,yt) ^ (a;o,yo)- 
By virtue of the equality 



(1) 
I — > 

(3) 



Xt-s 
Vt-s 



(2) 



Xt - Xt-s 

Vt - yt-s 



cos 9t sin 0t 
— sin 6t cos 9t 



)(;: 


- Xt-s 




' xl \ 




- yt-s 







reflection is a composition of the following transformations: (1) inversion of 

time on elastica; (2) reflection of the plane (x, y) in the center pc = {xt/2. yi/2) 
of the elastic chord I, i.e., the segment connecting its initial point (xq, yo) = (0, 0) 
and the endpoint {xt,yt), and (3) rotation by the angle {—Ot); see Fig. 27. 





Figure 27: Reflection of elastica in 
the center of chord Pc 



Figure 28: Reflection of elastica in 
the middle perpendicular 
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For reflection we have the decomposition 



(3) 



Vs 

COS 2x 



(1) 



Xt- 

Vt- 
sin 2x 



(2) 



Vt 



sin 2x cos 2% 

Xt - Xt-s \ (4) 

yt-s - Vt 



Vt 



+ 



COS Ut 

sin 9t 



— sm fft 
cos 9t 



— cos 2x — sin 2% 

— sin 2x cos 2x 

— cos 2x — sin 2x 

— sin 2x cos 2% 

Xt - Xt-s 

Vt-s - Vt 



)( 


Xt-s 




yt-s 


f Xt- 


:)] 


V Vt- 




xl \ 









where x is the polar angle of the point {xt,yt)- 

Xt 



cosx 



V^t + Vt 



sinx 



yt 



Vxi + 



yt 



Thus £^ acts on elasticae as a composition of 4 transformations: (1) inversion 
of time on elastica; (2) reflection of the plane {x, y) in the middle perpendicular 
l-^ to the elastic chord (3) and (4) rotations by the angles (— 2%) and 9t 
respectively; see Fig. 28. 

The symmetry acts on elasticae as reflection in the axis x. On the other 
hand, we have the following chain: 



(1) 



(2) 



cos 2x sin 2x 
sin 2x — cos 2x 

cos 2x sin 2x 
— sin 2x cos 2x 

,3 



COS 2x sin 2x 
sin 2x — cos 2x 



-ys 



ys 



this is a composition of: (1) reflection of the plane (x,y) in the elastic chord l\ 
and (2) rotation by the angle (— 2x); see Fig. 29. 



Figure 29: Reflection of elastica in 
chord / 



So, modulo inversion of time on elasticae and rotations of the plane {x,y), 
we have: 
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• is the reflection of elastica in the center of its chord; 

• is the reflection of elastica in the middle perpendicular to its chord; 

• is the reflection of elastica in its chord. 

9.6 Reflections of endpoints of extremal trajectories 

Now we can define the action of refiections in the state space M = K^ 
the action on endpoints of extremal trajectories: 



as follows: 



M^M, : qt^ ql, (9.9) 



£^ : I xc I I xtCos(9t + ytsin6it |, (9.10) 
-xt sinOt + yt cos Of 

Ot 

-2 " 




e : I Xt \ ^ I xtCos9t + yt sin 0t \ , (9-11) 
Xt sin 6t - yt cos 6t 



-3 ' 



: \ Xt \ ^ \ Xt I . (9.12) 

-yt 

These formulas follow directly from Propos. 9.1. Notice that the action of 
reflections £* : M — > M is well-deflned in the sense that the image £*(?() 
depends only on the point qt, but not on the whole trajectory {g^ | s e [0, t]}. 



9.7 Reflections as symmetries of the exponential mapping 

The action of reflections on the vertical subsystem of the normal Hamiltonian 
system (9.5) cicfincs the action of in the preimage of the exponential mapping 
by restriction to the initial instant s = 0: 

: u={l3,c,r)^v' = {l3\c\r), 

where {l3,c,r) — (/3o,co,r), {[3^,c^,r) — {jSQ.cly r) arc the initial points of the 
curves Vg = {Ps^Cs^r) and vl = {(3l,cl,r). The explicit formulas for (/3%c*) are 
derived from formulas (9.2)-(9.4): 

{(3\c') = {(3t,-ct), 
if3^c^) = {-f3t,ct), 
(/33,c3) = (-/3o,-co). 

So we have the action of reflections in the preimage of the exponential map- 
ping: 

e' : N^N, e\v) = v\ v, e N = T*^M. 
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Since the both actions of in N and M are induced by the action of £* on 
extremals As (9.7), we obtain the following statement. 

Proposition 9.2. Reflections e' are symmetries of the exponential mapping 
Expj : TV M, i.e., the following diagram is commutative: 



Expt 

N 



Expt Y 

N 



Expj 

iy\ -^qt 



Expj 



9.8 Action of reflections 

in the preimage of the exponential mapping 

In this subsection we describe the action of reflections 

in elliptic coordinates (Subsec. 7.3) in the preimage of the exponential map- 
ping N . 

Proposition 9.3. (1) If y = {k,(p,r) G A^i, then = {k,(fi^,r) £ Ni, and 

1 2K / 4K 
V +<^t = ^ mod^ 

2 ( ,4Ji:\ 
+ (yjt = I mod-^ I , 

, 2K ( AK 

ip — ip = —p= I mod 



(2) If V = (kjtjjjr) e N2, then = {k,ip^,r) G N2, moreover, 



and 



e N^, v"^ e N^, e N^, (9.13) 



V-^ + V-t = ^mod^^ 



tP'^ +tjj^ = Ol mod 



.2K 

( 2K 
- = mod 



(3) If v = {ip,r) e A^3, then = (<p%r) G N^, moreover. 
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and 

ip^ + ipt = 0, 
v'^ + <fit = 0, 

(^3 _ ^ = 0. 

(4) If V = {P,c,r) e Nq, then = {f3\d,r) S Nq, moreover, 
anrf 

(/3\ci) = (A,-c), 

(/32,c2) = (-/3„c), (9.15) 
(/33,c3) = (-/3,-c). 

Proof. Wc prove only item (1) since the other items axe proved similarly. 
The reflections preserve the domain A^i since 

: E E, s'^ : r r, 
e^, : c i— > — c, : c i— > c, 

this follows from equalities (9.2)-(9.3). Further, we obtain from (9.2) that 

= 0u = -ct, 

whence by virtue of the construction of elliptic coordinates (Subsec. 7.3) it 
follows that 

sii{s/r-^^) = sii{^ipt), cn(V?'</3^) = -cn{^/ript), 



thus (p^ + (ft ^ ^ (™°*^ V^j • "^^"^ expressions for action of the rest reflections 
in elliptic coordinates are obtained in a similar way. □ 

10 Maxwell strata 

10.1 Optimality of normal extremal trajectories 

Consider an analytic optimal control problem of the form: 

q = f{q,n), q&M, u € U. (10.1) 

q{0) = qo, q{h) = qu h flxed, (10.2) 

JtA(l,u]= <fi{q{t),u{t))dt ^min. (10.3) 
Jo 
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Here M and U are finite-dimensional analytic manifolds, and f(q,u), (p(q,u) 
are respectively an analytic vector field and a function depending on the control 
parameter u. Let 

hu{X) = {X,f{q,u))~^{q,u), X£T*M, q = tt{\) £ M, u&U 

be the normal Hamiltonian of Pontryagin Maximum Principle for this problem, 
see Subsec. 6.1 and [2]. Suppose that all normal extremals A* of the problem 
are regular, i.e., the strong Legendre condition is satisfied: 



hu{\t) <-S, 6> 0, (10.4) 

u{t) 



for the corresponding extremal control u{t). Then the maximized Hamiltonian 
H{\) = max/i„(A) is analytic, and there is defined the exponential mapping for 

time t: 

Exp, : AT = t; M ^ M, Exp, ( A) = tt o e*^ ( A) = . 

Suppose that the control u max;imizing the Hamiltonian /i„(A) is an analytic 
function u = u{\), A^G T*M. 

For covectors A, A G T*^M, we denote the corresponding extremal trajecto- 
ries as _ 
qs = Exp^(A), qs = Expg(A) 

and the extremal controls as 

u{s)=u(Xs), Xs = e'"(X). 

The time t Maxwell set in the preimage of the exponential mapping N = 
T*^M is defined as follows: 

MAX, = {a e TV I 3 A e TV : qs ^ qs, sG [0,t], qt = qt, Jt[q,u] = Jtiq,^}. 

J10.5) 

The inclusion A G MAX, means that two distinct extremal trajectories ^ qs 
with the same value of the cost functional Jt[q, u] = Jt[q, u\ intersect one another 
at the point qt ~ qt, see Fig. 30. 

The point qt is called a Maxwell point of the trajectory qs, s G [0,ti], and 
the instant t is called a Maxwell time. 

Maxwell set is closely related to optimality of extremal trajectories: such a 
trajectory cannot be optimal after a Maxwell point. The following statement is 
a modification of a similar proposition proved by S.Jacquet [15] in the context 
of sub-Riemannian problems. 

Proposition 10.1. If a normal extremal trajectory q^, s G [0,ti], admits a 
Maxwell point qt, t G (0,ti), then qs is not optimal in the problem (10.1)-(10.3). 
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Qt = qt, Jt[Q^u] = Jt[q,u\ 



Figure 30: Maxwell point qt 



Proof. By contradiction, assume that the trajectory qg, s G [0,ti], is optimal. 
Then the broken curve 

iqs, se[0,t], 
\qs, s e [t,ti] 



Is 

is an admissible trajectory of system (10.1) with the control 



^,^\u{s), sG[0,t], 
\u{s), s G [t,ti]. 

Moreover, the trajectory q'^ is optimal in the problem (10.1)-(10.3) since 

Jt^[q',u']= <fi{q'g,u'{s))ds= (p{q'g,u'{s)) ds + ifi{q'g,u' {s)) ds 
Jo Jo Jt 

= I ip{qs,u{s))ds+ I ip{qs,u{s))ds 
Jo Jt 

rti rti 

= Jt[q,u]+l ^{qs,u{s))ds = Jt[q,u]+ I if{qs,u{s))ds 
= Jtdl^i^ 



t Jt 



which is minimal since Qs is optimal. 

So the trajectory q'^ is extremal, in particular, it is analytic. Thus the 
analytic curves qg and q'^ coincide one with another at the segment s G [t,ti]. 
By the uniqueness theorem for analytic functions, these curves must coincide 
everywhere: qg = q'^, s £ [0,ti], thus q^ = q^, s £ [0,ti], which contradicts to 
definition of Maxwell point qt- □ 

Maxwell points were successfully applied for the study of optimality of 
geodesies in several sub-Riemannian problems [1,28]. We will apply this no- 
tion in order to obtain an upper bound on cut time, i.e., time where the normal 
extremals lose optimality, see [31-33] for a similar result for the nilpotent sub- 
Riemannian problem with the growth vector (2,3,5). 
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As noted in the book by V.I. Arnold [5], the term Maxwell point originates 
"in connection with the Maxwell rule of the van der Waals theory, according to 
which phase transition takes place at a value of the parameter for which two 
maxima of a certain smooth function are equal to each other" . 

10.2 Maxwell strata generated by reflections 

We return to Euler's elastic problem (3.4)-(3.9). It is easy to see that this prob- 
lem has form (10.1)-(10.3) and satisfies all assumptions stated in the previous 
subsection, so Propos. 10.1 holds for Euler's problem. 

Consider the action of reflections in the preimage of the exponential mapping: 

e' : N^N, e'{X) = X\ 

and denote the corresponding extremal trajectories 

g, = Exp,(A), g:=Exp,(A') 

and extremal controls (6.11) 

u{s) = Cs, u'{s) = c!,. 

The Maxwell strata corresponding to reflections e* are defined as follows: 

MAXj = {A e iV U: # qi = qt, Jt[q,u] = Jt[q\u']}, i = 1,2,3, t > 0. 

(10.6) 

It is obvious that 

MAXjcMAXt, i = 1,2,3. 

Remark. Along normal extremals we have 




In view of the expression for the action of refiections on trajectories of the 
pendulum (9.2)-(9.4), we have 

Jt[q\u'] = Jt[q,u], i = 1,2,3, 

i.e., the last condition in the definition of the Maxwell stratum MAX( is always 
satisfied. 

10.3 Extremal trajectories preserved by reflections 

In this subsection we describe the normal extremal trajcc;tories Qs such that 
ql = qs- This identity appears in the definition of Maxwell strata MAXJ (10.6). 

Proposition 10.2. (1) ql = qs <^ A^ = A. 
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(2) q'^ = qs ^ = A or A e N^. 

(3) = q, ^ A3 = A. 
Proof. First of all notice the chain 

qi = q, ei = es Pi-pl = Ps-Po, i= 1,2,3. (10.7) 

(1) Let ql = qs. By equality (9.2), (Sj = (it-s, thus we obtain from (10.7) 
that 

A_,-/?t = /3«-/3o. 
For s = t we have (3t = Po, thus 

Differentiating w.r.t. s and taking into account the equation of generalized 
pendulum (7.19), we obtain 

In view of equality (9.2), 

{(3lcl) = {(3,,cs) ^ {P\c') = {(3,c) ^ A=A\ 

Conversely, if A^ = A, then ql = qg. 

(2) Let ql = q^. In view of (9.3), /3f = -/?*-«, then (10.7) gives the identity 

-Pt-s + l3t = Ps-Po- 

Differentiating twice w.r.t. the equation of generalized pendulum (7.19), we 
obtain 

Ct-s = Cs => -rsm(3t-s = rsml3r (/3s = /Jo or = -/J^ or r = 0). 

If f3s = f3o, then Cg = 0, which means that X e N4U N^U Nr. If /3t-s = -/3s, 
then (/3s, Cs) = (/3s, Cs), thus A^ = A. Finally, the equality r = means that 
X€ NeUNr. So we proved that 

ql = q, => (A^ = A or A e uJ=47Vi). 

But if X e N4U N5U Nt, then /3s = or tt, Cs = (see Subsec. 8.1), and 
equality (9.3) implies that (/3g,c^) = {/3s, Cs), thus A^ = A. The implication 
=J» in item (2) follows. The reverse implication is checked directly. 

(3) Let q^ = qg. Equality (9.4) gives /3f = —/3s, and condition (10.7) implies 
that ps = /3o- Then Cs = 0. Consequently, A e N^U N^U N7. But if A G 

U A''5 U Nr, then A^ = A by the argument used above in the proof of item (2). 
The implication in item (3) follows. The reverse implication in item (3) 

is checked directly. □ 

Proposition 10.2 means that the identity ql = qg is satisfied in the following 

cases: 

(a) A' = A, the trivial case, or 

(b) A G iVg for i = 2. 
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10.4 Multiple points of the exponential mapping 

In this subsection wc study solutions to the equations g| = qt related to the 
Maxwell strata MAXj (10.6). 

Recall that in Subsec. 9.6 we defined the action of refiections in the state 

space M. Wc denote (j* = e^{q), q,q^ € M. 

The following functions are defined on M = ^ x up to sign: 



although their zero sets {P = 0} and {Q = 0} are well-defined. 
Proposition 10.3. (1) q^ = q 4^ 9 = (mod 27r), 

{2) q'^ = q <^ P = 0, 

{3)q^ = q <^ (y = and 61 = (mod tt)). 

Proof. We apply the formulas for action of reflections e' in M obtained in Sub- 
sec. 9.6. 

(1) Formula (9.10) means that 



If {x, y) = (0, 0), then q^ = {9, 0, 0) = g and P = 0, thus statement (2) follows. 
Suppose that {x,y) ^ (0,0), then we can introduce polar coordinates: 



^09 

P = xsm--ycos-, 

^ ■ 

Q = xcos- + j/sin-, 



: {9,x,y) 1-^ {—9,xcos9 + ysin9,—xs\ia9 + ycos9), 

which gives statement (1). 
(2) Formula (9.11) reads 

: q ~ {9, x,y) 1-^ q'^ = {9, X cos 9 + y sin 9, x sin 9 — y cos 9) . 



X = pcosx, 



y 



psmx. 



with p > 0. We have: 




cosYsin - — sin vcos - = <^ P = 0, 
^22 

and statement (2) is proved also in the case {x,y) ^ (0,0). 
(3) Formula (9.12) reads 

: q= {9, x, y) ^ q^ = {-0, x, -y), 
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thus 



y = -y 



e = (mod tt) 
y = 0. 



□ 



Notice the visual meaning of the conditions = qt for the corresponding 
arcs of Euler elasticae {xs,ys), s G [0,t] in the case x^+Pf ^ 0. As above, 
introduce the polar coordinates 

Xt = pt cos Xt, yt = Pt sin Xt 

with pt > 0. 

The condition 

qI = Qt et = o = eo <^ xt-do = xt- Ot 

means that the elastic arc has the same slope at its endpoints. The configuration 
corresponding to the inclusion qt G M^, where 



= {qeM\q^ =q} 



is shown at Fig. 31. 




■qt 




Figure 31: qt e Figure 32: qt e 

The condition 

qt =qt Pt=0 Xt-do = Ot- Xt 

means that the angle between the clastic arc and the elastic chord connecting 
{xo^Vo) to {xt,yt) reverses its sign. The configuration corresponding to the 
inclusion qt G M^, where 

M'^ = {qGM\q'^ = q} 



is shown at Fig. 32. 
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Finally, we have 



4 = Qt ^ 



?/ = and = (mod 2tt) 
y = and 9 = n (mod 27r) . 



Thus the set 



{qeM\q^ = q} 



has two connected components 



= {qe M \ y = 0, 6 = (mod 27r)}, 
M^- = {q M \ y = 0, 9 = 7: (mod 27r)}, 
= U M^- , U M^- = 0. 



See the illustrations to the inclusions qt e M^+, qt S at Figs. 33, 34 

respectively. 



It is easy to describe the global structure of the sets M\ The set = 
{q = {x, y,6) \ 9 = 0} is a. two-dimensional plane. It is the unique 2-dimensional 
Lie subgroup in the Lie group E(2) — the group of parallel translations of the 
two-dimensional plane Rx,y The set = {q = (a;, y,9) | P = 0} is the Mobius 
strip. Finally, and are straight lines. Notice that 

= M^n M^, 

M^- n = 0, M^- r\M'^ = {{9 = n,x = 0,y = 0)}. 

10.5 Fixed points of reflections 

in the preimage of the exponential mapping 

In order to describe fixed points of the reflections £* : N ^ N, we use elliptic 
coordinates {k, if, r) in A'' introduced in Subsec. 7.3. Moreover, the following 
coordinate will prove very convenient: 



While the values \/rip and \/r(pt correspond to the initial and terminal points 
of an elastic arc, their arithmetic mean t corresponds to the midpoint of the 
elastic arc. 



Figure 33: qt e 



Figure 34: qt e 



T = 



2 
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Proposition 10.4. Let v = {k,(fi,r) e Ni, then = £*(z/) = {k,(p'^,r) G Ni. 
Moreover: 

(1) i/i = J/ <^ cnr = 0, 

(2) v"^ = p <^ snr = 0, 

(3) = 1/ is impossible. 

Proof. We apply Propos. 9.3. The inclusion z/' e A^i holds. Further, 

1 1 2jr / , iK\ 

u = u <^ w = w <^ (yj + (Y3( = ^— mod —pr 
<^ T = K (mod 2ii:) <^ cnT = 0, 



= v 



V <^ = iy9 + (/?t = 0( mod 



4.K 



T = (mod 2;^") <^ snT = 0, 



<p3 = <^ 



2if 

7^ 



mod 1 which is impossible. 



□ 



Notice the visual meaning of the fixed points of the reflections £* : A^i — > A^i 
for the standard pendulum (7.1) in the cylinder (/3, c), and for the corresponding 
inflectional elasticae. 

The equality cnr = is equivalent to c = 0, these are inflection points of 
elasticae (zeros of their curvature c), see Fig. 35, 36. 









^^^^ ft 




Figure 35: cnr = 0, v G Ni 



Figure 36: Inflectional elastica cen- 
tered at inflection point 



The equality snr = is equivalent to /3 = 0, these are vertices of elasticae 
(extrema of their curvature c), see Fig. 37, 38. 
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Figure 37: snr = 0, u G Ni 



Figure 38: Inflectional elastica cen- 
tered at vertex 



In the domain A^2, we use the convenient coordinate 



T = 



2 



corresponding to the midpoint of a non-inflectional elastic arc. 

Proposition 10.5. Let v = {k,ijj,r) € N2, then z/' = £^{i^) = {k,4}\r) € N2- 
Moreover: 

(1) = 1/ is impossible, 

(2) v'^ = u <^ snrcnT = 0, 

(3) i/^ = 1/ is impossible. 

Proof. We apply Propos. 9.3. The inclusion i^^ € N2 holds. Implication (9.13) 
yields statements (1) and (3). We prove statement (2): 



Notice the visual meaning of the fixed points of the refiections e' : A^2 — > 
A^2- The equality snrcnr = is equivalent to the equalities /? = (mod tt), 
|c| = max;, min, these are vertices of non-inflectional elasticae (local extrema of 
their curvature c), see Figs. 39-42. 

Similarly to the previous cases, in the set A/3 we use the parameter 




<^ T = {mod K) <^ snrcnr=0. 



□ 



T = 



2 
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Proposition 10.6. Let v = {ip,r) € Ns, then = = {(p\r) e Ns. 

Moreover: 

(1) ly^ = v is impossible, 

{2) = v ^ T = 0, 

(3) = f is impossible. 

Proof. Follows exactly as in Propos. 10.5. □ 

The visual meaning of fixed points of reflection : ^ N^: the equality 
T = means that (3 = 0, \c\ = max, these are vertices of critical elasticae, see 
Figs. 43, 44. 




— TT c T 



Figure 43: t = 0, u e N3 Figure 44: Critical elastica centered 

at vertex 

Proposition 10.7. Let u = {(3,c,r) e N^, then = e\v) = {I3\c\r) e A^e- 
Moreover: 

(1) = u is impossible, 

(2) 1/2 = 1/ <^ 2/3 + ct = (mod 27r), 

(3) = u is impossible. 

Proof. Items (1), (3) follow from implication (9.14). Item (2) follows from (9.15) 
and the formula = Po + ct, see Subsec. 8.1. □ 

10.6 General description of the Mctxwell strata generated 
by reflections 

Now we summarize our computations of Maxwell strata corresponding to reflec- 
tions. 

Theorem 10.1. (1) Let v = [k,ip,r) € Ni. Then: 

I at = 9t I 6't = 
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(1.2) i/gMAX? <^ 

(1.3) z/eMAXj^ <^ 



v'^ J sn T 7^ 

ft =1t ^ {Pt = 

v'^^v ^ fy* = 

ql = qt 1 6»t = or TT. 



(2) Let V = {k, r) €N2. Then: 

(2.1) v e MAX(i ^ ^t = 



(2.2) V e MAX^ <^ 

(2.3) ^ 



Qi = Qt 

J sn r cn r ^ 
ql = qt ^ \Pt = 

ql = qt \0t = O orn. 



(3) Let u = {ip,r) e N3. Then: 

(3.1) u e MAXj^ <^ {^^^ ^ = 

[9t = Qt 

(3.2) .GMAX? ^ 1"'^'^ ^ 1"^' 

(3.3) . e MAX? ^ 1"'^'^ ^ 1^*=^ 

(4) MAXj r\Nj =9 for 1 = 1, 2, 3, j = 4, 5, 7. 
(6) Let veNe. Then: 

(6.1) 1/ e MAXj^ <^ S ! = 

(6.2) u e MAXj is impossible 

(6.3) . e MAX? o 1"'^'^ ^ 1^*=' 

* \ql = qt \9t = 0orn. 

Proof. In view of the remark after definftion of the Maxwell strata (10.6) and 
Propos. 10.2, we have 

MAXi = {ly G N \u' ^ i^, qi = qt}, i = l,3, 

MAX^ nNj = {u€ Nj \u^^u, q^ = qt}, j ^ 6, 

MAX^ DNq = 0. (10.8) 
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This proves the first implication in items (1.1) (3.3). The second implication 
in these items follows directly by combination of Propositions 10.4, 10.5, 10.6 
with Proposition 10.3. So items (l)-(3) follow. 

In the case v G 7V4 U A'^g U A^7 the corresponding extremal trajectory is 
{Xs,ys,6s) = (s,0,0), which is globally optimal since elastic energy of the 
straight segment is J = 0. By Propos. 10.1, there are no Maxwell points in 
this case. 

Finally, let v G Nq. Items (6.1) and (6.2) follow by combination of Pro- 
pos. 10.7 with Propos. 10.3. Item (6.2) was already obtained (10.8) from item 



Remark. Items (1.3), (2.3), (3.3.), (4), (6.3) of Th. 10.1 show that the Maxwell 
stratum MAX^ admits a decomposition into two disjoint subsets: 



In order to obtain a complete description of the Maxwell strata MAXJ, in 
the next section we solve the equations that determine these strata and appear 
in Th. 10.1. 

11 Complete description of Maxwell strata 

11.1 Roots of equation = 

In this subsection we solve the equation dt — that determines the Maxwell 
stratum MAX], see Th. 10.1. 



(2) of Propos. 10.2. 



□ 



MAX^ = MAX^+ UMAXf 



MAX^+nMAX^- = 0, 




We denote by 



A 
B 




, which denotes 



the condition A A B. 



Proposition 11.1. Let v = {k,(p,r) € A^i, then 



et = o ^ 



p = 2Kn, n e Z 
cnr = 0, 



where p = 



Vr{¥>t - f) _ y/rjift + y>) 
2 ' ^ 2 
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Proof. We have 

et=0 I3t = Po (mod 27r) <^ 



|dn(^^,,) = ±dn(v^(p) 



/3t /^o / , N 
— = — (mod tt) 

^ fsn(Vry)t) = sni^ip) 
[dn(^i^i) = dn(^i^) 

V^Vt = V^V (mod 4-^/4') 

y/ript = 1\prK — ^Jrip (mod i^^/rK) 



snp = 
cnr = 



p = 2Kn, n € 
cnr = 0. 



Proposition 11.2. Let v = {k,'il),r) G N2, then 

6t = <^ p= Kn, n e 



where p = 

z 

Proof. Let u € N:^, then 



et = o ^ 



(mod tt) 



2 2 

j sn(-v/rVt) = ± sn( Vr-'^) 

j^cn(-v/rVt) = ±cn(-v/ri/') 

V^^t = ^^ip (mod 4Y^/ir) 

= + (mod A^K) 

p = (mod 2ii:) 
p = /s: (mod 2K) 



p= Kn, n e 



If i/ G , then the same result is obtained by the inversion i : N2 — > N2 
Proposition 11.3. Let v € -/V3, then 

9^ = ^ t = 0. 

Proof Let u € N^, then 



= <^ 



(mod tt) <^ 



tanh(\/r(y3t) = ±ta,nh{^/r^p) 

^ =± ' 

cosh (y'r^t) cosh{^/r^p) 



□ 



□ 



{ cosh(v' 



sh{y/ript) cosh(v^(^) 
The same result is obtained for v e via the inversion i : A^3" — > 



t = 0. 



□ 
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Proposition 11.4. Let v £ Nq. Then 

et = <^ ct = 27rn, n e Z. 
Proof. We have 9t = ct in the case u G Nq. □ 

11.2 Roots of equation P = for u E Ni 

Using the coordinates 

r=^^^^ = V^(^+0, P=^^^^ = ^, (11.1) 
or, equivalent ly, 

= T - p, \frift = T + p, 

and addition formulas for Jacobi's functions (see Sec. 13), we obtain the follow- 
ing in the case v G N-^: 

sin — = 2fc sn p dn p en T ( dn^ t + cv? p sn^ r )/ , 

cos Y = (dn^P — k^sn^p cn^T)(dn^T + cv? p sn^r)/A^, 

xt = [2(2E(p) --p)(l - 2fc2snV) 

+ 8A;^ snp (cnp dnp (2E(p) — p) — snp + fc^ sn^p) cnr snr dnr 

-4fc2(p(-l+ sn2p(l-fc^(sn2p -2))) 

+ 2E(p)(l+ sn2p(-l + fc2(sn2p -2))) 

+ cnp snp dnp (— 1 + fc^ sn^ P (1 + cn^ r ))) sn^ r 

+ 2fc4(2 E(p) - p) sn2 p ( sn^ p - 2) sn^ r ]/(^/^A2) , 

t/t = 4fc[A;^cnp(2E(p) — p) sn^p cn^r snr dnr 

+ cnp dn^p snr ((p — 2E(p)) dnr + 2fc^ sn^p cnr snr) 

+ dnpsnp(2E(p)cnr(l + A;^(sn^p -2)sn^T) 

— pcnr(l + /c^(sn^p — 2)sn^T) 

+ dnr snr(l + fc^sn2p(sn2r - 2)))]/ {^/rA^), 

A = 1 — fc^ sn^ p sn^ r , 

^^^ 4fcsnrdnr A(p,fc) ^ . e iV„ (11.2) 
VrA 

/i(p, k) = snp dnp — (2 E(p) — p) cnp . 

In order to describe roots of the equation /i(p) = 0, we need the following 
statements. 

We denote by E{k) and K{k) the complete elliptic integrals of the first and 
second kinds respectively, see Sec. 13. 
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Proposition 11.5 (Lemma 2.1 [33]). The equatic 



2E{k) - K{k) = 0, 
has a unique root fcg G (0, 1). Moreover, 



fee [0,1), 



k e [0, fco) 
k e (A:o,l) 



2E-K >0, 
2E-K <Q. 



Notice that for fc — ^ we have 

V 2 



K = 



1 



r 



2K^+TT 



see [22], page 89, Chap. 3, exercise 24. Thus 

1 
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< fcn < 1. 



(11.4) 



The graph of the function k i-^ 2E{k) — K{k) is given at Fig. 45. Numerical 
simulations show that fcg ~ 0.909. 

To the value k = ko there corresponds the periodic Euler elastic in the form 
of figure 8, see Fig. 20. 

2E{k) - K(k) 

ir/2 




Figure 45: Definition of fco 



Proposition 11.6 (Propos. 2.1 [33]). For any k S [0, 1) the function 

/i(p, fc) = sn p dnp — (2 E(p) — p) cnp 
has a countable number of roots p\, n G Z. These roots are odd in n: 

P-n = -Pn: neZ, 
in particular, = 0. The roots p]^ are localized as follows: 
Pn e {-K + 2Kn, K + 2Kn), n E Z. 
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In particular, the roots p\ are monotone in n: 



Moreover, for n G N 



k e [0, ko) => p\ e {2Kn, K + 2Kn), 

k = ko ^ p\ ^ 2Kn, 
fce(fco,l) => p\e{-K + 2Kn,2Kn), 

where ko is the unique root of the equation 2E{k) — K{k) ~ 0, see Propos. 11.5. 

Proposition 11.7 (Cor. 2.1 [33]). The first positive rootp — p\ of the equation 
flip) — is localized as follows: 

k e [0,fcQ) ^ pie {2K,'SK), 

k^ko => p\ = 2K, 
fee (fco,l) ^ p} e {K,2K). 

The graph of the function k i— > p\{k) is shown at Fig. 46, and the graph k ^ 
p\{k) / K{k) is given at Fig. 47. Recall that the period of pendulum corresponds 
to p = 2K, this value is denoted at the axis of ordinates at Fig. 47. 

V\IK 




k- k 



0.6 . 



k 



Figure 46: k ^ p\, v € Ni 



Figure 47: k i-> p\/K, v e Ni 



Now we can obtain the following description of roots of the equation Pf = 
for V e Ni. 



Proposition 11.8. Let i/ G Ni. Then: 

Pf = ^ 
Proof. Apply Eq. (11.2) and Propos. 11.6. 



/i(p) = 
sn r = 



P = pi, neZ 
sn r =0. 



□ 
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Remark. We can propose a visual way for evaluating the roots to the equation 
fi{p) = 0, see Figs. 48, 49. Given an inflectional elastica, one should take its 
inflection point O, construct tangent lines to the elastica through the point O, 

and denote them AiBi, A2B2, ■ ■ ■ ; the tangents arc ordered by the length of the 
elastic arcs l(A„i?„). Then the number corresponds to the length 1(A„_B„); 

precisely, = — i„ = l(A„i?„), since elasticae are parametrized by arc 

length and in view of (11.1). 

On the arc AiBi the pendulum makes more than one oscillation in the case 
k < ko (Fig. 48), and less than one oscillation in the case k > ko (Fig. 49); 
thus in the first case pi > 2K{k), and in the second case pi < 2K{k). This 
observation provides one more illustration to Propos. 11.7 and Fig. 47. 
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) \ ' ) 
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Figure 48: Computing for k < ko Figure 49: Computing for k > ko 



11.3 Roots of equation P — for 1/ e N2 

In order to find the expression for Pt for u G N2 , we apply the transforma- 
tion of Jacobi's functions k ' ^ 
equality (11-2): 



Pt = 



^, see Subsubsec. 7.2.2 and (13.7), (13.8), to 

4fci sn(ri , fci ) dn(Ti , fci ) /i (pi , fci ) 



Vr(l - kf sn2(pi, /ci) sn2(Ti, /ci)) ' 



i/eATi, fcie(0,i), (11.5) 



Tl 



Pi 



2 



The both sides of equality (11.5) are analytic single- valued functions of the ellip- 
tic coordinates {ki, (p, r), so this equality is preserved after analytic continuation 
to the domain ki G (1,-|-(X)), i.e., v G N2 ■ 

Denote k2 = then the formulas for the transformation fc 1— > ^ of Jacobi's 
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functions (13.7), (13.8) give the following: 

p ^ 4^^'Hn.^)cu(n.^).A(ju.^) 

* ^^^(l-^sn2(pi,i)sn2(n,i)) 

44fc2sn(g^,fc2)cn(g^,fc2)/2(P2,fe2) 
V^(l-^fcisn2(|i,fc2)fcisn2(fi,fc2)) 
4sn(T2, fca) cn(T2, fc2)/2(p2, ^2) 



where 



Vr(l - /e| sn2(p2,A;2)sn2(T2,fc2))' 



= fc2 = 2fc2 = 2 ' ("-^^ 

Pi %/^(y't - y) \/^(^'/ - c^) , „x 



and 

h{P2,k2) = fi [ pi, 



1 

sn ( p„ 1) dn (p„ ^) - (2E (p„ 1) - P:) cn (p„ 1 
fc2sn ( ^.^2] cn f ^.^'2 



1 



= — [^2 Sn(p2,fc2) Cn(p2,fc2) +dn(p2,fc2)((2 - /C2)j32 - 2E(p2, fe))]. 
K2 



Summing up, we have 

_ 4snr cnr/2(p, fc) 
* ~ Pa ' 



e iV2+, (11-8) 



/2(P,fc) = ^[fc'snpcnp + dnp((2-fc2)p-2E(p))], 

\A^(V't + V') VriA-tp) A 1 ;2 2 2 
r = ^ , p = , A = l — fcsnpsnr. 

We will need the following statement. 
Proposition 11.9 (Proposition 3.1 [33]). The function f2{p) has no rootsp ^ 0. 
Proposition 11.10. Let v € N2. Then 



<^ 



snr CUT =0. 
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Proof. In the case u G N2, we obtain from (11.8) and Propos. 11.9: 

Pt = ^ 



/2(P) = 

sn T cn r =0 



p = 

snr cnr =0. 



The case v G N2 is obtained by the inversion i : N2 — > A^2 • The inversion 
i maps as follows: 

(/3,c,r) 1-^ -c,r), {0,x,y) ^ {-9,x, -y), 



thus equality (11.8) yields 



4 sn T cn T /2 (p, k) 



and the statement for the case v G Ny follows. 



(11.9) 

□ 



11.4 Roots of equation P = for v E 

Passing to the limit k 1—0 in equalities (11.8), (11.9), we obtain the following: 

4tanhT/2(i5, 1) 



Pt = ±- 



•yrcoshT(l — tanh^ptanh^ r) 



Hp A) = , lim /2(p,fc) 

K— >1 — 



Vr{v>t + y) 

2 

2p — tanhp 



coshp 



(11.10) 



(11.11) 



Proposition 11.11. Let v G N^. Then 

ft = ^ 



J5 = 

r = 0. 



Proof. We have {2p — tanhp)' = 2 — ^^^j^2 ^ > 1, thus 

/2(p, 1)=0 <^ 2p-tanhp = <^ P = 0, 
and the statement follows from (11.10). □ 

11.5 Roots of equation P = for v E Nq 
Proposition 11.12. If v € Nq, then Pt = 0. 

Proof. Pt = Xt sin — — Vt cos — = - sm sin — 1 1 — cos ct) cos — = 0. □ 

■' 2^ 2c 2c^ ^2 

The visual meaning of this proposition is simple: an arc of a circle has the 
same angles with its chord at the initial and terminal points. 
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11.6 Roots of system y — 0, 6 — 

Notice that 



0t = O 

yt = o 



Ot = 

Pt = xt sin ^ -yt cos ■ 



0, 



so wc can replace the first system by the second one and use our previous results 

for equations 6t = and Pt = 0. 

Proposition 11.13. Let u € Ni. Then 



et = o 

Pt = 



k = ka, p = 2Kn, 

P = Pn, CUT =0, 

p = 2Kn, snr =0, 



n e z. 



Proof. By virtue of Propos. 11.1, 11.8, we have 



9t = 
Pt = 



I p = 2Km or cn r =0 
\p = p]^ov sn T =0 



I p = 2Km 

\p = pi 



or 



I cnr =0 

\p = pi 



or 



I p = 2Km 
snr =0 



or 



(11.12) 

cnr = 
snr =0. 
(11.13) 



By Propos. 11.6, 



I p = 2Km 

\p = pi 



P = Pn = ^Kn <^ 



\k = ko 
I p = 2Kn. 



Now it remains to notice that the system cnr =0, snr = is incompatible, 
and the proof is complete. □ 



Proposition 11.14. Let u G N2. Then 



et = o 

Pt=0 



p = Kn, T = Km, 
p = 0, 



n, m &Z. 



Proof. Taking into account Propos. 11.2 and 11.10, we obtain 



Pt=0 



I p — Kn 

I p = or r = Km 



p = Kn, T = Km, or 
p = 0. 



□ 



Proposition 11.15. Let v G N^- Then 

\Pt = Q 



t = 0. 
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Proof. Follows immediately from Propos. 11.3, 11.11. □ 
Proposition 11.16. Let v e Then 

= 
= 

Proof. Follows immediately from Propos. 11.4, 11.12. □ 
11.7 Roots of system y — 0, 9 — n for u e Ni 

The structure of solutions to the system yt = 0, dt = tt is much more complicated 
than that of the system = 0, 0^ = studied above. 
First of all, for any normal extremal 



et = o 



cos — = 

^ e, 9, (11-14) 

Qt = xt cos — + t/t sm — = 0. 



From now on we suppose in this subsection that v & N\. 
In the same way as at the beginning of Subsec. 11.2, in the coordinates r, p 
given by (11.1) we obtain 

cos ^ = ( dn^ p — sn^ p cn^ r ) ( dn^ t + cn^ p sn^ t ) / 

= (l-2fc^sn^p + fc^sn^p sn^T)(dn^T + fc^ cn^p sn^ r )/A^, 
Qt = 2E(p) -p + fc^sn^r(2cnpsnpdnp - (2E(p) -p)(2- sn^p)). 



Thus 



cos^=0 <^ sn^r = (2fc^sn^p - l)/(fc^sn2T). 



Substituting this value for sn^ r into Qt, we get rid of the variable r in the 
second equation in (11.14): 

2 

'3tlsn2T=(2/c2sn2p-l)/(fc2sn2^) = gU^ p ' 

gi{p,k) = {1 — k'^ + k^ cn^ p){2E{p) — p) + cnp snp dnp{2k^ sn^ p — 1). 

(11.15) 

So we can continue chain (11.14) as follows: 

cos^=0 ^ fsn^r = (2fc2sn2p -l)/(fc2sn2T) 
Qt = \5i(p,fc) = 0. 

We proved the following statement. 
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Proposition 11.17. Let v e Ni. Then 

f sn2 T = (2fc2 sn2 p - l)/{k'^ sn^ r ) 




. , (11.16) 
gi{p,k) = 0. 

Now wc study solvability of the second system in (11.16) and describe its 
solutions in the domain {p G {0,2K)}. For the study of global optimality of 
normal extremal trajectories, it is essential to know the first Maxwell point. By 
Propos. 11.1, the first Maxwell point corresponding to occurs at p = 2K, so 
for the study of the global optimal control problem we can restrict ourselves by 
the domain {p G (0, ^K)}. What concerns the local problem, in the forthcoming 
paper [34] wc show that only the Maxwell strata MAX], MAX^ , but not MAXf 
are important for local optimality. But for the global problem, the stratum 
MAX( is very important: in fact, on this stratum extremal trajectories lose 
global optimality, i.e., MAXj provides a part of the cut locus [34]. 

The second system in (11.16) is compatible iff the equation gi{p, k) = has 

, ■ i ,^ , , ^ 2fc^ sn^;> - 1 ^ , . , . 

solutions [p, k) such that < — — ^ < 1, or, which is equivalent, 

k'^ sn'^ p 

2A;2sn2p - 1 > 0. (11-17) 
After the change of variable 



r dt 

p = F{u,k)= ^ <^ u = am{p,k), 

Jo VI - fc2 sin2 t 



(11.18) 



where am(p. A;) is Jacobi's amplitude (see Sec. 13), we obtain 

gi{p, k) = hi{u, k), 

hi{u, k) = {l-k^ + fc2 cos'' u){2E{u, k) - F{u, k)) 

+ cosusinu\/ 1 — k'^ sin^ u{2k^ sin^ u — 1). (11.19) 

Denote 

hi{u, k) 



h2{u, k) 



1 — + fc2 cos^ u 

cosusinw\/l — fc2 sin^ u(2A;^ sin^ u — 1) 



= 2E{u, k) - F{u, k) + 



k'^ + k'^ cos^ u 



(11.20) 

ai{u,k), (11-21) 



a direct computation gives 

d /i2 sin^ u\/2 — fc^ + fc^ cos 2u 
du ~ 4\^(l-fc2 + A;2cos4u)2 

ai{u, k) = Co + ci cos2u + C2 cos^ 2u, (11.22) 
Co = 8 - lO/c^ + 4A:^ 
ci =4fc2(3-2fc2), 
C2 = 2fc2(2A:2 - 1). 
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First we study roots of the function ai(u, k). In view of the symmetry relations 

ai(u + TT, fc) = ai{iT — u,k) = ai{u, k), (11.23) 

we can restrict ourselves by the segment u G [0, |]. Consider the corresponding 
quadratic function 



ai = cq + cit + C2t^ , t = cos2u & [—1,1]. 



li k = then 



ai=A{l + t)=0 <^ t = -l. 



(11.24) 



Let k e 1). Then cq > 0, ci > 0, C2 > 0, thus ai > for t e [0, 1]. On 
the other hand, 

ailt=_i = Co - ci + C2 = 8(1 - A;^)(l - 2k^) < 0. 

Thus the quadratic function ai = cq + cit + C2t^ has a unique root ta^ at the 
interval t G (—1,0). Consequently, the function ai{u,k) given by (11.22) has a 
unique zero = | arccosf^^ at the segment u G [0, |], moreover, € ( |, f ). 
We prove the following statement. 

Proposition 11.18. (1) The set ^^{u,k) 
smooth curve. 



e JR X 



ai{u, k) = o| *s a 



(2) There is a function 



such that 



V2' 



'Wai(fc) 



2' 



A; e 



72' 



and for k = —j= , 1 
v2 



ai(u, A;) = 44> u = Ua-^{k) + nn = — + nn, 



while for k € 



Moreover, 



V2' 

ai{u,k) = <^ 



U = Mai (^) + 2"''^ 

U = TT — Uai {k) + 2-Kn. 



72" 



72' 



(11.25) 



(11.26) 
(11.27) 
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Proof. We assume in this proof that (u, A;) e K x 
(1) We have 



— — = —2 sin 2u(ci + 2c2 cos 2u) . 
ou 

It is easy to show that 

{u,k)\ai{u,k) =Q, ^(m,A;)=0 

= l^fc) = (1 + '^"':^) ' = (1 + '^"'-^ 

Further, for u = ^ + im we have 

ai = Co - ci + C2 = 8(1 - k'^){l - 2fc2), 
which has regular zeros at fc = = 1. Thus at the points (w, fc) = 

+ TTO, and {u,k)= (f + 7rn, l) we have 

ai{u,k) = 0, ^^{u,k)^0. 

By imphcit function theorem, the equation ai {u, k) = determines a smooth 
curve. 

(2) For /c G ( '-'-)' already defined before this proposition 



1 / TT 

Uai{k) = -arccostai(fc) e (^-, 



where tai{k) G (—1,0) is the unique root of the quadratic polynomial ai = 
cq + cit + C2t^. We define now 

For k = we have by virtue of (11.24): 

TT 

ai = t = cos 2u = — 1 -s^- w = — + Tzn. 

For A; e (757 1) and m G [O, |] , we get 

Oi = ^ t = COS2m = e ( — 1, 0) ^ U=Uai, 

and in view of the symmetry relations (11.23), implication (11.26) follows. 
Let k = 1, then 

ai(w, 1) = 2 + 4cos2w + 2cos^ 2u = <^ «=^+7rn, 

and implication (11.25) is proved. 

Finally, the regularity relations for the function Uai{k) specified in (11.27) 
follow from the implicit function theorem by the argument of item (1). □ 
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The plot of the curve {ai{u, k) = 0} in the domain {u E (0, tt)} is presented 
at Fig. 50. 




Figure 50: ai{u, k) = 



The distribution of signs of the function oi (u, k) in the connected compo- 
nents of the domain {ai{u, k) ^ 0} shown at Fig. 50 follows from the relations 



u = 



ai(u, fc) = Co + ci + C2 = 8 > 0, 



ai{u, fc) = 0, u =^ — + 27rn 



dai 
du 



^0. 



Now we study the structure and location of the curve 

' 1 



7/11 = i {u,k) e (0,7r) X 



= ,1 



I hi{u,k)^0}, 



so below in this subsection we suppose that u G (0,7r), k e 1]. 
Recall that the function 

h2{u, k) = hi{u, k) {l-k^ + k^ cos^ u)~^ 

> 

has derivative (11.21) 

d /i2 sin^ uV2^^1?2~+~P~cos2u 



du 4^/2(1 -fc2 + fc2cos4y)2 



ai(M, fc). 



>o 



A direct computation from (11.19) gives 

2 



/ii(u, fc) + o(m^), m ^ 0, 
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thus 



hi{u, k) > and h2{u, k) > as u ^ +0. 
If M G (0, (A:)), then 

dh2 



(11.28) 



ai{u, k) > 



du 
hi > 0, 



>0 



/i2 T w.r.t. u 



h2>0 



thus 7/1, n{ue (0,Uai(fc))} = 0. 

Now we study the curve 7/1, in the domain {u G [uai{k), |]}. We have 



hi{u, fc) = (1 - k^){2E{k) - F(fc)), 
hiiu,k) = 2{2E{k) - F{k)). 



(11.29) 



Notice that the unique root fco of the equation 2E(k) — F{k) = satisfies 
the inequality fco G 1), see (11-4). 

Taking into account Propos. 11.5, we obtain: 



k e 



or M = TT 



hi{u,k) > 0, 
k = ko, (^u = ^ 01 u = TT^ hi{u, k) = 0, 

k e (fcoi 1) , = or u = TT^ /ii(u, fc) < 0. 

If fc € (fco, 1) then: 

U = Uai{k) /i2(w, fc)>0, 

TT 



U G 



Woi(fc), 



/i2(u, fc) J, W.r.t. u, 

/l2(w,fc) < 0, 



, TT - ?iai (fc) 



=> /i2(w, fc) i W.r.t. u, 

U = 'IT — Uai{k) /l2(w, fc)<0, 

u G [tt — Uai(fc),7r] /12(m, fc)T W.r.t. u, 

U = 'K /l2(w, fc)<0. 

Consequently, for fc G (fco, 1) the equation /i2(w, fc) = 0, or, equivalently, hi{u, fc) = 
0, has a unique root u = ^/(^(fc) at the interval u G (0,7r), moreover, Uh2{k) G 

K(fc),f). 

A similar argument shows that for fc = fco the equation /i2(w, fc) = has a 
unique root u = Ufi^{ko) at the interval u G (0,7r), moreover, ^/(^(fco) = f . 
In particular, we proved that 



fc G [fco, 1) 



h2{n -Ua,{k),k) <0, /ii(7r - (fc), fc) < 0. (11.30) 
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Now we determine the largest root of the function 

a{k) = hiin - Uai{k),k), fce(^,l). 



a/2' 

Notice that implication (11.30) means that a{k) < for G [ko, 1). By virtue 

thus the continuous function a{k) has roots at the interval k G (^^, A'o j . Denote 

fc* = supj/c e ^-^,A;o^ I a(A;) =o| , (11.31) 

see Fig. 51. It follows that fc* € (^^i^o^ 
If fc e (fc,,l), then: 

u = Uai{k) hi{u,k)>0, 
u = Tr — Uai{k) => hi{u,k)<0. 

Thus there exists a function 

Ufei : (fc*,l) -> , u = UhAk), 

such that for k £ (fc*, 1), u £ (0, tt — ^^^(fc)) 

{u,k)€jhi ^ hi{u,k)^0 u = Uhi{k), 

(fc) < Uhi (fc) < TT - Wai {k)- 

Further, we define 

= TT - Mai(fc*), Uhi{k*)=U*, (11.32) 

M?ii(l) = Mai(l) = |- 

For k G (/c*, 1), u e (0, tt), the curve 'jhi does not intersect the curve {ai{u, k) = 
0}. Taking into account equalities (11.20), (11.21), we conclude from implicit 
function theorem that u^^ e C[/c*, 1] n C°°(A;,, 1). 
We proved the following statement. 

Lemma 11.1. There exist numbers fc, e (;^) ^o)> e (f , ^) anrf a function 

Uh, ■■ [fc*,l] ^ (^1'^ 

SMc/t that for k € [fc,, 1],- u G (0, tt — (fc)) 

/ii(u, fc)=0 <^ u = Uhi{k). 
Moreover, uu^ e C[fc,, 1] n C°°{k», 1) and Uhi(fc*) = u*, Uhi{ko) = u{l) = f . 
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u 



k 



Figure 51: Definition of fc« 



Figure 52: Plot of fc i— > Uh^{k) 



A plot of the function Uhi{k) is presented at Fig. 52. 
Numerical simulations give fc* k, 0.841, k, 1.954. 

Lemma 11.1 describes the first solutions in u to the equation hi{u,k) ~ for 
fee [A:*, 1]. Now we show that this equation has no solutions for k e l:^^ j^*); 
u G (0,7r]. First we prove the following statement. 

Lemma 11.2. The function h2{u,k) defined by (11.20) satisfies the inequality 



dk 128fc(l-/c2)(l_fc2_^fc2cos4u)2' 

Ci = C2-E(u, k) + c^Fiu^ k) + C4^di sin 2u, 
C2 = -2{2k'^ - 1)(8 - 5fc2 + fc2(4cos2u + cos4u))2 < 0, 
C3 = -2(1 - fc2)(8 - 5fc2 + /c2(4cos2'u + cos4u))2 < 0, 
C4 = V4 - 2fc2 + 2fc2 cos 2m > 0, 

di = 72 - 90fc2 + 28^;-* - (48 - 97k^ + 34fc^) cos 2m + (8 - Qk^ + 4/0"*) cos 4m 
+ fc^(2fc^ - 1) cos 6m, 

so it is enough to prove that di{u^ fc) > for k £ (-4^, 1), m € , We have 



ei = 16 - 21fc2 + Qk'^ _ (12 - 25fc2 + 10fc'')t + (4 - S/c^ + 2k^)t^ + fc2(2fc2 _ 1)^3^ 
t = cos 2m e [-1,0], 

so it suffices to prove that for any k € 1); ^^'^ cubic polynomial ei(t) is 
positive at the segment i G [—1,0]. At the boundary of this segment we have: 




Proof. Direct computation gives 
dh2 ci 



di = 4ei, 



t = ei 16-21fc2 + 6fc'* > for fc2 e [0,1], 

t=-l ei = 16(2 + fc2)(l + fc2) >0, 
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thus it is enough to prove that 



te{-l,0), -^=0 ei>0. (11.33) 



We have 



: -12 + 25fc2 - lOfc^ + (8 - 6A;2 + Ak'^)t + (6A;^ - 3fc2)^^ 

:0 <^ t = ti{k) or t = t2{k), 



dei 

_ -4 + 3A;2 - - 2^4 - ISP + 43A;4 - 48/c6 + 16A;8 
^ ~ 3A;2(2P-1) ' 

-4 + 3fc2 - 2fc^ + 2\/4 - 15A;2 + 43/;;4 - mfi + 16fc8 

i2(A;) = 



3A;2(2fc2 _ 1) 

27k^{2k^ - 1)2 

ai{k) = (1 - fc2)(8 - 37fc2 + 146/c'' - 250fc^ + 224A;^ - 64fc^°), 
/?i(fc) = 4 - 15fc2 + 43^" - 48fc^ + 16A:^ 

and it is enough to prove that a\ > /Sf for k £ (:^) !)• 
We have 

as = a? - (3f = 27fc''(l - /c2)2(l - 2k^fa4, 
ai = 7 - 19fc2 + 71A;4 - 32A;^ 

Via standard calculus arguments, one can easily prove that for k e {^A) 

wc have a4 > 0, thus as > and ai > Pi^'^. The statement of this lemma 
follows. □ 

Let fcG i^--^)-- (f'T")- implicit function theorem, we obtain from 
Lemma 11.2 that the equation h2{u,k) = defines a smooth curve such that 
each its connected component is the graph of a smooth function k = k{u), 
^^e(f,f). 

We have /i2(it*,/c») = 0, so there exists such a connected component con- 
taining the point (u*, fc*); denote by A; = kh2 (u), u € (|, ^) the function whose 
graph is this component: 

h2{u,kh2{u)) = 0, kh^iu*) = K- 

Notice that 

[TT 

k=kh2{u) <;4> u = Uhi{k), ke[k^,ko],ue 
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Now wc prove that there are no other connected components of the curve 
{/i2(w, k) = 0} in addition to the component given by the graph of k = kh2{u). 
By contradiction, suppose that there is such a component k = kh^iu), u e 
(f > T")- "^^^ curve {k = kh2{u)} cannot intersect the curve {u — 'K—Ua^{k) = 0} 
at a point where k > k^, (this would contradict definition (11.31) of the number 
A;*) or = A;* (this would mean that kh^iu) = kh2{u))- Consequently, 



TT 37r 
2'T 



(11.34) 



Consider the limit lim„^^/2+o kh2 {u) ■ 

a) Suppose that there exists k = lim„_>^/2+o ^/i2(^)- By virtue of inequal- 
ity (11.34), k e [^,ko]. Since h2{u,kh2{u)) = 0, it follows that /i2(f ,fc) = 0, 
thus /ii(f ,fc) = 0. By virtue of (11.29), it follows that 2E{k) - F(k) = 0. By 
Propos. 11.5, we have k = ko, which is impossible since A^ijd) = ko and the 
curve {h2{u, k) = 0} is smooth at the point (fc = /cq, u = f ). 

b) Consequently, the limit lim„_>^/2+o (^) does not exist. But we can 
choose a converging sequence (u„, kn) — > (f + 0, k), k G [:^,kQ], and come to a 
contradiction in the same way as in item a). 

So we proved that the curve 



TT 37r 
2'T 



1 



X I -^,fco ) I h2{u,k) = 



{u, k) e 

consists of the unique connected component 

k = kh^iu), u G 

We have ^-^^ = and in view of equality (11.21) and Lemma 11.2, 

d kh 

it follows that sgn ^ = sgnai{u,k). Thus 

du 



TT 37r 
2' T 



u e 



T 



ai{u,kh2{u)) < 

ai{u, kh^iu)) > 



dkh2_ 
du 

dkh. 



<0, 

^ >o, 



so fc = fc* is the minimum of the function k = kh^ {u) , u G ( ^ , 
Then it follows that 



4 J 



(u, k) e [0 

Taking into account that 



37r 



1 
7^ 



h2{u,k) > 0. 



(u, k) e 



37r 



^/2 



ai{u, k) > 



dh2 
du 



>0, 
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we obtain finally 

(u, k) e (0, tt] X 



V2 



/i2(w, fc)>0 ^ A:)>0. 



In particular, the pair is the unique solution to the system hi{u,k) = 

ai{u, fc) = in the domain fc) € (|, ^] x 1^ • 
Summing up, we proved the following statement. 

Proposition 11.19. (1) The set = {{u, k) € (0, tt] x (0, 1] | hi{u, k) = 0} 
is a smooth connected curve. 

(2) The system of equations hi{u,k) = 0, ai{u,k) = has a unique solution 
(w*,fc*) G (f , ^,1 • Moreover, fc* G [^^ka^, u., G (f , ^). 

(3) The curve jhi does not intersect the domain {(w, fc) G (0,7r] x [0, A:*)}. 

(4) There exist functions 

u = Uh,{k), Uft, e C[fc*,l] nC°°(fc*,i), 
k = kh,{u), kh,eC° 

such that in the domain {(u, fc) G (0,7r] x [fc*, 1]} holds the following: 
^hi n {m e (0,M*], k G [fc*, 1]} = {u = Uh^k)}, 



2" 



,fc e 



V2 



7/,! n |w G 

T/ie function u = Uh-^ (fc) satisfies the hounds: 
k G [fc*,fco) Uh^k) G 

TT 

fc = fco UhAk) = ^, 



{fc = kh^u)}. 



2'T 



fcG(fco,l) ^ G (^,|) . 

In particular, for fc G [fc*, 1] 

min{u > | h\{u,k) = 0} = Uhi{k). 

Now we return to the full system (11. IG), in particular, to the condition of 
compatibility (11.17). After the change of variable (11.18) this condition reads 

Piu,k) = 2fc2sin2w - 1 > 0. 

We prove that this inequality holds on the curve u = Uh^k), fc G [fc*, 1]. 



84 



We have: 



(3{u, k) ^ sin^ u 



1 

2P 



^ cos 2u = 



fc2 



^ ai(u, = — — ^ > for fc e [fc*, 11. 

In other words, the curve {/3(u, k) = 0} is contained in the domain {ai{u, k) > 
0}. Thus the function /3(u, k) preserves sign on each connected component of 
the domain Da-^ = {ai(u, fc) < 0, fee [fc*, 1]}; since fc) = 2k'^ — 1 > 0, the 
function /3(m, k) is positive on Da^ . On the other hand, the curve {u — Uh^ (fc) | 
A; € [fc*, 1]} is contained in the domain Da^, thus 



k e [/c*, 1], w = Mhi(fc) 



/3(u,fc) > 0. 



(11.35) 



The plot of the curves {hx = 0}, {ai = 0}, {/3 = 0} is presented at Fig. 53. 
The elastica corresponding to fc = fc* is plotted at Fig. 54: for this elastica, the 
tangent hne at the inflection point touches the preceding and the next waves of 
the elastica; moreover, fc = A;* is the minimal of such k. 




y/2 



/ii = 



Figure 53: The curves {hi = 0}, {ai = 0}, {/? = 0} 



We return back from the variables (u, k) to the initial variables (p, k) via the 
formulas (11.18), and obtain the following statement. 

Proposition 11.20. Let the function gi{p,k) be given by (11.15). 

(1) The set 

jg, = {(p, k)\ke{OA), pe (0, 2K{k)), gi{p, k) = 0} 

is a smooth connected curve. 
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Figure 54: Elastica with k = 



(2) The curve does not intersect the domain 

{{p,k)\k€iO,k.), p€{0,2K{k))}. 

(3) The function 

P = Pg,{k)^F{Ul{k),k), Pg^€C°°{k,,l), 

satisfies the condition 

mm{p > I gi{p,k) = 0} =Pg^{k), k e [A;,,l). 
The function p = Pg^ (fc) satisfies the hounds: 

k e [k^,ko) ^ PgAk) G (^K, ^K^ , 
k = ko Pgi {k) = K, 

fce(fco,i) PgAk)&(^lK,Ky 

(4) For any k e [k^,l) 

p = PgM 2fc2sn2(p,fc)-le(0,l]. 

(5) If k € (0, fc*), then the system of equations (11.16) has no solutions (p, r) 
with p G {0,2K{k)). If k £ [A;,,l), then the minimal p G {0,2K{k)) such 
that the system (11.16) has a solution {p,t) is p = Pg^{k). 

So we described the first solution to system (11.16) derived in Propos. 11.17. 
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11.8 Roots of system y — 0, 9 — n for u e N2 

Similarly to Propos. 11.17, we have the following statement. 
Proposition 11.21. 

(a { 2 2sn2p-l 

r*:" ^ k-sn^p (11.36) 

where 

g\ip,k) = y[fc^cnp snp dnp{2sv? p — 1) 

+ (1 - 2sn2p + fc2sn4p)(2E(p) - (2 - fc^)^)]. (11.37) 

Proof. Let S A^2~- We apply the equivalence relation (11.14). Further, in 
order to obtain expressions for cos ^ and Qt through the variables T2, P2 given 
by (11.6), (11.7), we apply the transformation of Jacobi's functions A; i— > ^ in 
the same way as we did in Subsec. 11.3, and obtain 

6t [1 -2sn^ p + k'^sn^ p sn^ t){c!1^ T + dv? p sn^ t) , „„x 

cos - = ^2 ' (11-38) 

Qt = ^[(2E(p)-(2-A2)p) 

+ snV(2fc^cnp snp dnp - (2E(j5) - (2 - fc^)p)(2 - fc^sn^p)]. 

Consequently, cos ^ = Q -S^- sn^ r = ^Tg,^"^ • Then direct computation 
gives 

2 

sn2T=(2sn2p-l)/(fe2sii2p) = ^2 g^2 p 9l {Pi k) , 

where the function gi{p,k) is defined in (11.37). The statement of this propo- 
sition is proved for v G iV^, and for £ iV^ it is obtained via inversion 
i : N+ ^ N^. □ 

Now we study solvability of the system of equations (11.36) in the domain 
p e {0,K). This bound on p is given by the minimal p = K hi points in MAX^, 
see Propos. 11.2. 

After the change of variables (11.18), we have 

gi{p,k) = hi{u,k) 
_ 1 
~ k 



k^ cos u sin u\/ 1 — fc^ sin^ u{2 sin^ u — 1) 
+ (1 - 2 sin^ u + k^ sin* u){2E{u, k) - {2 - k'^)F{u, k)) 

pe{o,K) ^ "e(o,|). 
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Introduce the functions 



/i2(w, k) = -777— A;) = 1 — 2sin^ u + k"^ sin^ u. 

P{u,k) 



Then 

dhn P 



:sin^(2w)/i3(u,fc), (11.39) 



9u 16/3 Vl - k^sin^u 
hsiu, k) = 2(2 -k^ + 2k^ cos 2u + (2 - k^) cos^ 2m). 

The quadratic polynomial 

/i3(t) = 2(2 - fc^ + 2A;^i + (2 - /c^)i^) 

is positive for all t e M., thus h^{u,k) > for all u G M, fc e (0,1). Then 
equality (11.39) iniphes that > 0, thus /12 increases in u on all intervals 
where /?(«,, /e) 7^ 0. It is easy to show that 

/3(u, A;) = 0, u e ^0, u = ■U/3(fc) = arcsin --^^^=, k' = -\/l — fc^. 

If w e (0, up{k)), then 

/3(w, A;)>0 =J> /i2(u, fc)t in u 
and since /i2(0,A;) = hi{k) = 0, 

/i2(«,fc)>0 /ii(M,fc)>0. (11.40) 

If M = Ui3{k), then 



/^iK^)= (I^^>0- (11-41) 



Further, 



u^U0ik) + O =^ /3(w,fc) ^ -0, /ii(m,A;) ^ ^^-^^ >0 
=J> ft.2(w, k) —00. 
Finally, if u e (■u^(fc), |] , then 

P{u,k)<0 => h2{u,k)'\ in u. 

We have 

'vr ,\ 1 



/..(|,fc)=^7(fc), 
7(A:) = 2.E(A;) - (2 - k'^)F{k). 
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By Propos. 11.23 presented below and proved in [33], 7(fc) < for any k G (0, 1). 
Thus 

/i2(|,fc)<0 =J> h2{u,k)<0 y u € (u(}{k),'^ 

hi{u,k)>0 V u e (u/3(fc), I . (11.42) 
Summing up inequalities (11.40), (11.41), (11.42), we proved that 
hi{u,k)>o yue(o,^ , fee (0,1). 

We return back to the initial variable p via the change of variables (11.18), 
and obtain the following statement. 

Proposition 11.22. Let the function gi{p,k) be given by (11.37). Then for 

any k G (0, 1), p £ (0, K{k)) we have gi{p, k) > 0. 

In fact, numerical simulations show that the equation gi{p, k) = has solu- 
tions p > K. 

Here we present the statement used above in the proof of Propos. 11.22. 

Proposition 11.23 (Lemma 2.3 [33]). The function jik) = 2E{k) - (2 - 
k'^)K{k) is negative for k G (0, 1). 

11.9 Roots of system y = 0, 9 = n for f E 

Proposition 11.24. If v £ N3, then the system of equations yt = 0, Ot — n is 
incompatible for t > 0. 

Proof Let v e A^3+. We pass to the limit fc ^ 1 - in Propos. 11.21, 11.22 

and obtain that the system of equations yt = 0, = tt has no roots for p G 
{0,K{l-0)),p= ButK{l-0)=\imk^i_oK{k)=+oo. Thus the system 
in question is incompatible for t > and v G N^. The same result for 1/ G 
follows via the inversion i : — > . □ 

11.10 Roots of system y = 0, 9 = n for u G Nq 

Proposition 11.25. If v G Nq, then the system of equations yt = 0, 9t = is 
incompatible. 

Proof. As always, we can restrict ourselves by the case u e Nq. Then it is 
obvious that the system is incompatible: 

1 — cos ct 

yt = =0, Vt = ct = n + znk. 

c 

□ 
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11.11 Complete description of McLxwell strata 

Now we can summarize our previous results and obtain the following statement. 

Theorem 11.1. (1.1) Ni n MAX^^ = {ueNi\p= 2Kn, cnr ^ 0}, 

(1.2) Ni n MAXj^ = {uGNi\p = pi, snr ^ 0}, 

(1.3+) 7VinMAXf+ = {ueNi\ {k,p) = {ko,2Kn) or {p = p^.cnr = 0) or {p = 
2Kn, SILT = 0)}, 

(1.3-) NiDMAX^- = {i^eNi\gi{p,k) = 0, sn^ t = {2k^ sn^ p -1) / {k^ sn^ p)y , 
iVinMAX3-n{p e (0,2K)} = {k G [fc,,l), p = pg,{k), sn^r = 
(2A:2sn2p - l)/(fc2 sn^p)}^ 

(2.1) N2 n MAXt^ = {veN2\p = Kn, cnrsnr 7^ 0}, 

(2.2) N2 n MAX^ = 0. 

(2.3+) iVa n MAX-^+ = {i^ e N2 \ p = Kn, snr cnr = 0}, 

(2.3-) N2nMAX^- = {ly e N2 \gi{p,k) = 0, sn^ r = (2sn2p -l)/(A;2sn2p)}3, 

N2 n MAX'^- n{p e (0, if)} = 0, 

(3.1) 7V3 n MAXt^ = 0, 

(3.2) N3 n MAX^ = 0. 
(3.3+) N3 n MAXf+ = 0, 
(3.3-) N3 n MAXf- = 0, 

(6.1) Ne n MAXj^ = {u e Nq \ ct = 2-Kn}, 

(6.2) A/-6 n MAXj^ = 0, 

(6.3+) Nq n MAXf+ = {v &N(,\ct = 2-Kn}, 

(6.3-) Nq n MAXj^- = 0. 

Proof. It remains to compile the corresponding items of Th. 10.1 with appro- 
priate propositions of this section: 

(1.1) Propos. 11.1, 

(1.2) Propos. 11.8, 

(1.3+) Propos. 11.13, 

(1.3-) Propos. 11.17, 11.20, 

■'where the function gi(p, k) is given by (11.15) 
^where fc» and pg-^(k) arc described in Propos. 11.20 
^where the function gi(p,k) is given by (11.37) 
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(2.1 
(2.2 
(2.3+ 
(2.3- 
(3.1 
(3.2 
(3.3+ 
(3.3- 
(6.1 
(6.2 
(6.3+ 
(6.3- 



Propos. 11.2, 
Propos. 11.10, 
Propos. 11.14, 
Propos. 11.21, 11.22, 
Propos. 11.3, 
Propos. 11.11, 
Propos. 11.15, 
Propos. 11.24, 
Propos. 11.4, 
item (6.2) of Th. 10.1, 
Propos. 11.16, 
Propos. 11.25. 

□ 



12 Upper bound on cut time 

Let Qs, s > 0, be an extremal trajectory of an optimal control problem of the 
form (5.9)-(5.11). The cut time for the trajectory is defined as follows: 

tcut = sup{ti > I is optimal on [0,ti]}. 

For normal extremal trajectories = Expjj(A), the cut time becomes a function 
of the initial covector A: 

tcut : N = T;^M ^ [0, +oo], t = icut(A). 

Short arcs of regular extremal trajectories are optimal, thus tcut (A) > for 
any X G N. On the other hand, some extremal trajectories can be optimal on 
an arbitrarily long segment [0, ti], ti € (0, +oo); in this case tcut = +oo. 

Denote the first Maxwell time as follows: 

tY^^{\) = inf{t > I A e MAXJ. 

By Propos. 10.1, a normal extremal trajectory Qg cannot be optimal after a 
Maxwell point, thus 

tcut(A)<tr^(A). (12.1) 

Now we return to Euler's elastic problem. For this problem we can define 
the first instant in the Maxwell sets MAX\ i = 1,2, 3: 

t^^^\X) = inf{t > I A e MAXj}. 
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Since tf^^(A) < tf^^'(A), we obtain from inequality (12.1): 
tcut(A) <min(tfA^*(A)), i=l,2,3. 

Now we combine this inequality with the results of Sec. 11 and obtain an upper 
bound on cut time in Euler's elastic problem. To this end we define the following 
function: 

t : N ^ (0,+oo], Ah^t(A), 
2 

A e iVi ^ t = ^pi(fc), 

™i„w.).ri,*))^{^;;];_'' 11^^ (12^2) 

2k 

AeTVa ^ t=^pi{k), p,{k) = Kik), 

XeNe t = — r, 

|c| 

A e -/Vs U iV4 U A/5 u A/'t t = +00. 

Theorem 12.1. Let Xg N. We have 

icut(A) < t(A) (12.3) 

in the following cases: 

(1) A = {k,p, t) G A^i, cnr snr 7^ 0, or 

(2) X&N\Ni. 

Proof. (1) Let A = {k,p,T) G iVi, cnr snr 9^ 0. Then Th. 11.1 yields the 
following: 

k e (0, fco] ^ t(A) = ^2if = 

Vr 

fce(fco,l) ^ t(A) = 4W(fc) = tr'''(A). 
(2) Let A = {k,p,T) e iV2, then we obtain from Th. 11.1: 

snr cnr ^0 ^ t(A) = = ^max^ (A), 

snr cnr =0 ^ t(A) = = fr""'^ (A). 

Vr 

If A = (/3,c,r) e A'e, then Th. 11.1 implies that 

t{X) = Pr=tf^''\X)=tf^''\X). 
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If A e TVa, then there is nothing to prove since t(A) = +00. 
If A e A'4 U -/V5 U A^7, then there is also nothing to prove since in this case the 
extremal trajectory is optimal on the whole ray s G [0, +00), and tcut(A) = 



In the proof of Th.12.1, we used the explicit description (12.2) of the function 
Pi{k) = niin(2if (fc),p}(A:)) which follows directly from Prop. 11.5. 

Remark. In the subsequent work [34] we prove that if A = {k,p,T) G A^i and 
cnr snr =0, then the corresponding point qt = Expj(A), t = t{X) is conjugate, 
thus the trajectory qs is not optimal for s > t(A); consequently icut(A) < t(A), 
compare with item (1) of Th. 12.1. So the bound (12.3) is valid for all A G VV. 

Notice the different role of the Maxwell strata MAX^+ and MAX^~ for 
the upper bound of the cut time obtained in Th. 12.1. On the one hand, the 
stratum MAX^"*" generically does not give better bound on cut time than the 
strata MAX^, MAX^ since generically 



see Th. 11.1, this follows mainly from the fact that the system of equations de- 
termining the stratum MAX"^^ consists of the equations determining the strata 



see Th. 10.1. 

The situation with the stratum MAX^~ is drastically different. By item 
(1.3-) of Th. 11.1, we have 



t(A) = +00. 



□ 




MAX^ and MAX^: 




u = 



{k,p,T) eiViHMAX^ 



(12.4) 



fcG[fc*,l), p = pg^{k), sn^T = 



2k'^ sn^p 



- 1 



G [0,1]. 



(12.5) 



fc^ sn^ p 



Moreover, from Propos. 11.20, Th. 11.1, Propos. 11.6 it follows that 



fce[fc*,fco) pg,{k) < -K <2K = pi{k), 

k = ko ^ pg,{k) = K <2K = pi{k), 
k€{kQ,l) ^ pg,{k) <K <pl{k)=pi{k). 



That is, Pgi{k) < pi{k) for all k G [fc*, 1), consequently. 



(A) < t(A) = min (ff (A), (A)) 



for all A = z/ defined by (12.4), (12.5). 

It is natural to conjecture that for such A we have 




(A) 



(12.6) 
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and wc will prove this equality in the subsequent work [34] . 

Although, the covectors X = v defined by (12.4), (12.5) form a codimension 2 
subset of N, so the equality (12.6) defines the cut time for a codimension one 
subset of extremal trajectories. The question on exact description of cut time 
foe arbitrary extremal trajectories is under investigation now. 

An essential progress in the description of the cut time was achieved via the 
study of the global properties of the exponential mapping. Moreover, a precise 
description of locally optimal extremal trajectories (i.e., stable Euler elasticae) 
was obtained due to the detailed study of conjugate points. These results will 
be presented in the subsequent work [34]. 

13 Appendix: 

Jacobi's elliptic integrals and functions 

We base upon the textbooks of D.F. Lawden [22] and E.T. Whittaker, G.N. Wat- 
son [38]. 

13.1 Jacobi's elliptic integrals 

Elliptic integrals of the first kind: 




and of the second kind: 




Complete elliptic integrals: 






Jo 



13.2 Definition of Jacobi's elliptic functions 



<^ = amu ^ u = F{ip,k) 



dnu= \/l — k"^ sn^ u, 
E(u) = E{amu,k). 



cnu = cosamzt. 



snu = smamw, 



(13.1) 
(13.2) 

(13.3) 
(13.4) 
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13.3 Standeird formulas on Jacobi's elliptic functions 
Derivatives with respect to u 

am' u = dnu, 
sn' u = cnudnu, 
cn' u = — sn ?i dn u, 
dn' u = —k^ snucnu, 

Derivatives with respect to modulus k 

d snu 1 , k 2 1 T^/ \ 1 

— -- — = —ucnudnu+ -rrsnucn u——— -rr^ E(u) cnudnu, 

ok k 1 — fc^ fc(l — k'') 

d cnu 1 , k 2 1 T^/ ^ , 

— — — = —-U cnu dnu — — r sn ucnu + tt, snu dnu, 

ok k 1 — k-' k{l — k^) 

d dnu ^ 2 1 , T^/ N 

— — — = r^sn udnu — Kusnucnu H — Jifu) snucnu, 

ok 1 — k'' 1 — k-^ 

9 E(u) k , ,2 ^ T^/ X 2 

— — — = — 7 snucnudnu — Kusn u — r tfu) cn u, 

ok 1 — k^ 1 — fc^ 

dK _ E-{1- k^)K 

Ilk ~ k{l - fc2) ' 

dE E-K 



dk k 
Integrals 

Addition formulcis 

sn(u + v) 

cn(u + v) 
dn(u + v) 



/ dn^ tdt = E{u). 
Jo 



sn u cn w dn u + cn u dn u sn u 
1 — fc^ sn^ u sn^ V ' 
cn u cn u — sn u dn u sn u dn u 



1 — fc^ sn^ u sn^ V ' 
dn u dn w — /c^ sn u cn u sn u cn u 
1 — k^ sn^ usn^ v ' 
E(u + u) = E(u) + E(u) — fc^ sn u sn u sn(u + v) . 



Degenerat ion 

A; ^ +0 snu^sinu, cnu^cosu, dnu^l, E(u) u, 

(13.5) 

fc ^ 1 — => snu^tanhu, cnu, dnu ^ — \ — , E(u)-^tanhu. 

cosnu 

(13.6) 
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Transformation fc 

k 




/ 1\ 1 /M ,\ 1 - fc^ 
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